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ABSTRACT 



Radiation images of an object constituted of tissues 
exhibiting different levels of radiation absorptivity are 
formed with two kinds of radiation having different 
energy levels. A first original image signal SOI and a 
second original image signal S02 representing the radi- 
ation images are detected. A superposition image signal 
is generated with a weighted addition expressed as 

S~ (N2 2 /(N1 2 + A3 2 )) X SOI + (W 2 /(M 2 +N2 2 ))- 
KSOZ 



wherein Nl and N2 respectively represent noise con- 
tained in the first and second original image signals. A 
first image signal representing a first image primarily 
composed of patterns of first tissues of the object is 
generated by subtracting the first and second original 
image signals from each other. A first smoothed image 
signal representing a first smoothed image, in which 
noise of the first image has been reduced or eliminated, 
is generated by smoothing the first image signal. A 
second image signal representing a second image pri- 
marily composed of patterns of second tissues of the 
object is generated by subtracting the superposition 
image signal and the first smoothed image signal from 
each other. 

56 Claims, 11 Drawing Sheets 
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METHOD FOR FORMING ENERGY 
SUBTRACTION IMAGES 

BACKGROUND OF THE INVENTION 
L Field of the Invention 

This invention relates to a method for forming an 
energy subtraction image wherein, from a plurality of 
radiation images, an energy subtraction image is formed 
which includes little noise and which has good image 
quality and can serve as an effective tool in, particu- 
larly, the efficient and accurate diagnosis of an illness. 
This invention also relates to a superposition processing 
method and apparatus for radiation images, wherein an 
addition process is carried out on two image signals 
representing radiation images of a single object 

2. Description of the Prior Art 

Techniques for reading out a recorded radiation 
image in order to obtain an image signal, carrying out 



10 



15 



even when the energy intensity of the radiation to 
which the stimulable phosphor sheet is exposed varies 
over a wide range. More specifically, since the amount 
of light which the stimulable phosphor sheet emits 
when being stimulated varies over a wide range and is 
proportional to the amount of energy stored thereon 
during its exposure to the radiation, it is possible to 
obtain an image having a desirable density regardless of 
the energy intensity of the radiation to which the stimu- 
lable phosphor sheet was exposed. In order for the 
desired image density to be obtained, an appropriate 
read-ont gain is set when the emitted light is being de- 
tected and converted into an electric signal (image sig- 
nal) to be used in the reproduction of a visible image on 
a recording material, such as photographic film, or on a 
display device, such as a CRT. 

In the radiation image recording and reproducing 
systems wherein recording media, such as X-ray film or 
stimulable phosphor sheets, are used, subtraction pro- 



appropriate image processing on the image signal, and 20 cessing techniques for radiation images are often carried 



then reproducing a visible image by use of the pro- 
cessed image signal have heretofore been known in 
various fields. For example, as disclosed in Japanese 
Patent Publication No. 61(1986>5 193, an X-ray image is 
recorded on an X-ray film having a small gamma value 25 
chosen according to the type of image processing to be 
carried out, the X-ray image is read out from the X-ray 
film and converted into an electric signal, and the elec- 
tric signal (image signal) is processed and then used for 
reproducing the X-ray image as a visible image on a 30 
photocopy, or the like. In this manner, a visible image 
having good image quality with high contrast, high 
sharpness, high graininess, or the like, can be repro- 
duced. 

Also, when certain kinds of phosphors are exposed to 35 
radiation such as X-rays, a-rays, £-rays, y-rays, cathode 
rays or ultraviolet rays, they store part of the energy of 
the radiation. Then, when the phosphor which has been 
exposed to the radiation is exposed to stimulating rays, 
such as visible light, light is emitted by the phosphor in 40 
proportion to the amount of energy stored thereon 
during its exposure to the radiation. A phosphor exhibit- 
ing such properties is referred to as a stimulable phos- 
phor. 

As disclosed in US. Pat. Nos. 4,258,264, 4,276,473, 45 
4,315,318, 4,387,428, and Japanese Unexamined Patent 
Publication No. 56(1981>11395, it has been proposed to 
use stimulable phosphors in radiation image recording 
and reproducing systems. Specifically, a radiation 
image of an object, such as a human body, is recorded 
on a sheet provided with a layer of the stimulable phos- 
phor (hereinafter referred to as a stimulable phosphor 
sheet). The stimulable phosphor sheet, on which the 
radiation image has been stored, is then scanned with 
stimulating rays, such as a laser beam, which cause it to 
emit light in proportion to the amount of energy stored 
during exposure to the radiation. The light emitted by 
the stimulable phosphor sheet, upon stimulation thereof, 
is photoelectrically detected and converted into an 
electric image signal. The image signal is then used 60 
during the reproduction of the radiation image of the 
object as a visible image on a recording material such as 
photographic film, on a display device such as a cath- 
ode ray tube (CRT), or the like. 

Radiation image recording and reproducing systems 65 
which use stimulable phosphor sheets are advantageous 
over conventional radiography using silver halide pho- 
tographic materials, in that images can be recorded 



50 



55 



out on image signals detected from a plurality of radia- 
tion images of an object which have been recorded on 
the recording media. 

With the subtraction processing techniques for radia- 
tion images, an image is obtained which corresponds to 
a difference between a plurality of radiation images of 
an object recorded under different conditions. Specifi- 
cally, a plurality of the radiation images recorded under 
different conditions are read out at predetermined sam- 
pling intervals, and a plurality of image signals thus 
detected are converted into digital image signals which 
represent the radiation images. The image signal com- 
ponents of the digital image signals which represent the 
image information recorded at corresponding sampling 
points in the radiation images are then subtracted from 
each other. A difference signal is thereby obtained 
which represents the image of a specific structure or 
part of the object represented by the radiation images. 

Basically, subtraction processing is carried out with 
either the so-called temporal (time difference) subtrac- 
tion processing method or the so-called energy subtrac- 
tion processing method. In the former method, in order 
for the image of a specific structure (for example, a 
blood vessel) of an object to be extracted from the 
image of the whole object, the image signal represent- 
ing a radiation image obtained without injection of 
contrast media is subtracted from the image signal rep- 
resenting a radiation image in which the image of the 
specific structure (for example, a blood vessel) of the 
object is enhanced by the injection of contrast media. In 
the latter method, such characteristics are utilized that a 
specific structure of an object exhibits different levels of 
radiation absorptivity with respect to radiation with 
different energy levels. Specifically, an object is ex- 
posed to several kinds of radiation with different energy 
levels, and a plurality of radiation images are thereby 
obtained in which different images of a specific struc- 
ture are embedded. Thereafter, the image signals repre- 
senting the plurality of the radiation images are 
weighted appropriately and subjected to a subtraction 
process in order to extract the image of the specific 
structure. The applicant proposed novel energy sub- 
traction processing methods using stimulable phosphor 
sheets in, for example, U.S. Pat Nos. 4,855,598 and 
4,896,037. 

A plurality of radiation images, which are subjected 
to energy subtraction processing, will herein be referred 
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to as the "original images". An image signal represent- 
ing a subtraction image is obtained by subtracting the 
image signals representing the original images from 
each other. Therefore, the image signal representing the 
subtraction image has a lower signal-to-noise ratio (S/N 
ratio) than the image signals representing the original 
images. As a result, the problems occur in that the 
image quality of the subtraction image becomes worse 
than the image quality of the original images. 

By way of example, energy subtraction processing is 
often carried out in the manner described below. Specif- 
ically, an object, such as the chest of a human body, 
which is constituted of soft tissues and bones, is exposed 
to several kinds of radiation with different energy lev- 
els, and a plurality of radiation images of the object are 
thereby obtained. The plurality of the radiation images 
are read out, and a plurality of image signals represent- 
ing the radiation images are generated. Energy subtrac- 
tion processing is then carried out on the plurality of the 
image signals. From the energy subtraction processing, 20 
a soft tissue image signal is obtained which represents a 
soft tissue image primarily composed of patterns of the 
soft tissues of the object. Alternatively, a bone image 
signal is obtained which represents a bone image pri- 
marily composed of patterns of the bones of the object. 
Thereafter, the soft tissue image is reproduced as a 
visible image from the soft tissue image signal, or the 
bone image is reproduced as a visible image from the 
bone image signal. In the soft tissue image, the patterns 
of the bones have been erased. Therefore, patterns, 
which were behind the bone patterns or were rendered 
imperceptible by the bone patterns in the original im- 
ages, become more perceptible in the soft tissue image 
than in the original images. Also, in the bone image, the 
patterns of the soft tissues have been erased. Therefore, 
patterns, which were behind the soft tissue patterns or 
were rendered imperceptible by the soft tissue patterns 
in the original images, become more perceptible in the 
bone image than in the original images. Accordingly, a 
subtraction image can be obtained which is well 40 
matched to the purposes of diagnosis. However, be- 
cause the soft tissue image and the bone image are ob- 
tained from the subtraction processing, the problems 
occur in that noise components have been emphasized 
in the soft tissue image and the bone image than in the 45 
original images. From this point of view, the image 
quality of the soft tissue image and the bone image 
could not heretofore been kept good. 

Accordingly, in pending U.S. patent application Ser. 
No. 654,450, the applicant proposed a method for form- 50 
ing an energy subtraction image wherein a subtraction 
image is formed in which noise has been reduced. The 
proposed method comprises the steps of: 

i) after a plurality of radiation images of an object are 
recorded on recording media by irradiating several 55 
kinds of radiation with different energy levels to 
the object, which is constituted of a plurality of 
tissues exhibiting different levels of radiation ab- 
sorptivity with respect to the several kinds of radia- 
tion with different energy levels, and a plurality of 60 
original image signals representing the plurality of 
the radiation images are then detected, 

generating a first image signal, which represents a 
first image primarily composed of patterns of first 
tissues of the object, from the plurality of the origi- 65 
nal image signals, 

ii) generating a first smoothed image signal by pro- 
cessing the first image signal, the first smoothed 
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image signal representing a first smoothed image in 
which noise components of the first image have 
been reduced or eliminated, and 
iii) generating a second image signal by subtracting 
the first smoothed image signal from an original 
image signal, the second image signal representing 
a second image primarily composed of patterns of 
second tissues of the object. 
In pending U.S. patent application Ser. No. 654,450, 
the applicant also proposed a method for forming an 
energy subtraction image wherein, after the aforesaid 
first process is carried out in order to generate the first 
image signal, the second process is carried out in order 
to generate the second image signal. Thereafter, a third 
process is carried out, which comprises the steps of: 

a) generating a second smoothed image signal by 
processing the second image signal, the second 
smoothed image signal representing a second 
smoothed image in which noise components of the 
second image have been reduced, and 

b) generating a new first image signal by subtracting 
the second smoothed image signal from an original 
image signal, the new first image signal represent- 
ing a new first image primarily composed of the 
patterns of the first tissues of the object. 

The proposed method for forming an energy subtrac- 
tion image may be modified such that the second pro- 
cess or the third process is repeated, and noise compo- 
nents of the image signal may be reduced even further. 

With the methods for forming an energy subtraction 
image, which are proposed in pending U.S. patent appli- 
cation Ser. No. 654,450, the image signal components of 
the first image signal and the second image signal are 
added to each other which represent the image informa- 
tion stored at corresponding picture elements in the two 
radiation images. In this manner, a superposition image 
signal SO is obtained, which can be expressed as 

SO={SO\+SOl)/2 

wherein SOI represents the first image signal, and S02 
represents the second image signal. In the superposition 
image signal SO, noise components of the image signal 
have been reduced. Therefore, the addition process is 
advantageous for the subsequent processes. 

However, the superposition image signal is obtained 
with the formula shown above wherein the values of 
the first image signal and the second image signal are 
merely averaged. Therefore, the degree of reduction in 
the noise components is not necessarily be appropriate. 
Thus a need exists for a method for forming an energy 
subtraction image, which yields an energy subtraction 
image in which the noise components have been re- 
duced even further. 

SUMMARY OF THE INVENTION 

The primary object of the present invention is to 
provide a method for forming an energy subtraction 
image, wherein an energy subtraction image is formed 
in which noise has been reduced to a higher extent and 
which has better image quality and can serve as a more 
effective tool in, particularly, the efficient and accurate 
diagnosis of an illness than an energy subtraction image 
obtained with a conventional method for forming an 
energy subtraction image. 

Another object of the present invention is to provide 
a superposition processing method for radiation images, 
wherein weight factors appropriate for two sets of 
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image signals are employed during superposition pro- 
cessing carried out on the two sets of the image signals, 
and a superposition image having good image quality is 
thereby obtained from the superposition processing. 

The specific object of the present invention is to 5 
provide an apparatus for carrying out the superposition 
processing method for radiation images. 

The present invention provides a first method for 
forming an energy subtraction image comprising the 
steps of: 10 

i) after two kinds of radiation images of an object are 
formed with two kinds of radiation having differ- 
ent energy levels, the object being constituted of a 
plurality of tissues exhibiting different levels of 
radiation absorptivity with respect to the two kinds 15 
of radiation having different energy levels, 

detecting first and second original image signals rep- 
resenting the two kinds of the radiation images, 

ii) generating a superposition image signal by carry- 
ing out a weighted addition expressed as 20 

S= (A^VCM 2 +/V2 2 )) +(MV<M 2 +/V2 2 ))- 
XS02 

wherein SOI represents the first original image 
signal, Nl represents noise contained in the first 25 
original image signal, S02 represents the second 
original image signal, and N2 represents noise con- 
tained in the second Original image signal, 

iii) generating a first image signal, which represents a 
first image primarily composed of patterns of first 30 
tissues of the object, by subtracting the first origi- 
nal image signal and the second original image 
signal from each other, 

iv) generating a first smoothed image signal by 
smoothing the first image signal, the first smoothed 35 
image signal representing a first smoothed image in 
which noise components of the first image have 
been reduced or eliminated, and 

v) generating a second image signal by subtracting 
the superposition image signal and the first 40 
smoothed image signal from each other, the second 
image signal representing a second image primarily 
composed of patterns of second tissues of the ob- 
ject 

In the above-described first method and also in the 45 
below-described various other methods according to 
the present invention including those defined in claims, 
the expression of "with two kinds of radiation having 
different energy levels" does not necessarily means two 
separate radiations but include two kinds of radiations 50 
originated from a single radiation wherein, for example, 
one is a direct radiation from an X-ray source and the 
other is a radiation from the same source passing 
through a phosphor sheet and/or a filter or the like the 
low energy components of which are filtered out. 55 
Therefore, the two kinds of radiation images can be 
formed one after another using different radiations of 
different energy levels or formed simultaneously using a 
single radiation with two recording media or stimulable 
phosphor sheets placed one upon another with or with- 60 
out a filter interposed therebetween. When the filter is 
not used, the stimulable phosphor placed one another 
functions as a filter to filter out the low energy compo- 
nents. 

The first method for forming an energy subtraction 65 
image in accordance with the present invention may be 
embodied in various, substantially identical manners. 
For example, each of the steps of the first method for 



forming an energy subtraction image in accordance 
with the present invention may be divided even further 
into a plurality of steps. Alternatively, the operations 
may be carried out in different orders. 

By way of example, in cases where the object is a 
human body, the term "a plurality of tissues exhibiting 
different levels of radiation absorptivity" as used herein 
means various tissues of the human body, such as bones 
and soft tissues, or mammary cancer tissues and mam- 
mary gland tissues. In cases where the object is an inor- 
ganic body, such as a machine subjected to non-destruc- 
tive inspections, the term "a plurality of tissues exhibit- 
ing different levels of radiation absorptivity*' as used 
herein means portions of the inorganic body, which are 
constituted of different materials, and the like. 

The present invention also provides a second method 
for forming an energy subtraction image comprising the 
steps of: 

i) after two kinds of radiation images of an object are 
formed with two kinds of radiation having differ- 
ent energy levels, the object being constituted of a 
plurality of tissues exhibiting different levels of 
radiation absorptivity with respect to the two kinds 
of radiation having different energy levels, 

detecting first and second original image signals rep- 
resenting the two kinds of the radiation images, 

ii) generating a superposition image signal by carry- 
ing out a weighted addition expressed as 

xsoz 

wherein SOI represents the first original image 
signal, Nl represents noise contained in the first 
original image signal, S02 represents the second 
original image signal, and N2 represents noise con- 
tained in the second original image signal, 

iii) carrying out a first process for generating a first 
image signal, which represents a first image primar- 
ily composed of patterns of first tissues of the ob- 
ject, by subtracting the first original image signal 
and the second original image signal from each 
other, 

iv) thereafter carrying out a second process, which 
comprises the steps of: 

a) generating a first smoothed image signal by 
smoothing the first image signal, the first 
smoothed image signal representing a first 
smoothed image in which noise components of 
the first image have been reduced, and 

b) generating a second image signal by subtracting 
the superposition image signal and the first 
smoothed image signal from each other, the sec- 
ond image signal representing a second image 
primarily composed of patterns of second tissues 
of the object, and 

v) thereafter carrying out a third process, which 
comprises the steps of: 

a) generating a second smoothed image signal by 
smoothing the second image signal, the second 
smoothed image signal representing a second 
smoothed image in which noise components of 
the second image have been reduced, and 

b) generating a new first image signal by subtract- 
ing the superposition image signal and the sec- 
ond smoothed image signal from each other, the 
new first image signal representing a new first 
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image primarily composed of the patterns of the of the first to fourth methods for forming an energy 

first tissues of the object subtraction image in accordance with the present inven- 

An image having better image quality can be obtained tion is modified such that the superposition image signal 

by repeating the second and third processes in the sec- i s generated by carrying out a weighted addition ex- 

ond method for forming an energy subtraction image in 5 pressed as 

accordance with the present invention. Specifically, the 

present invention further provides a third method for SMSttS(Ski+skZ))xSOt+(Skx/(Skt+sia)y 

forming an energy subtraction image comprising the xSOl 
steps of: 

i) carrying out the processes in the second method for 10 wherein SOI represents the first original image signal, 
forming an energy subtraction image in accordance Ski represents a sensitivity of the first original image 
with the present invention, and signal, S02 represents the second original image signal, 

ii) thereafter repeating the following once or several and Sk2 represents a sensitivity of the second original 
times: image signal. 

a) a new second process for generating a new sec- 15 By way of example, in cases where an original image 
ond image signal by carrying out the second signal is detected from a stimulable phosphor sheet, the 
process in which the new first image signal ob- term "sensitivity of an original image signal" as used 
tained from the third process is taken as the first herein means the sensitivity of the original image signal, 
image signal in the second process, the new sec- In cases where an original image signal is detected from 
ond image signal generated by the new second 20 X-ray film, on which an X-ray image has been recorded, 
process representing a new second image pri- the term "serisitivity of an origind image signal" as used 
marily composed of the patterns of the second herein means the sensitivity of the X-ray film, 
tissues of the object, and In the methods for forming an energy subtraction 

b) a new third process for generating a new first image in accordance with the present invention, the 
image signal by carrying out the third process in 25 superposition image signal is generated from the 
which the new second image signal is taken as weighted addition expressed as 

the second image signal in the third process, the 

new first image signal generated by the new 5= (A3 2 /(/vi 2 + A2 2 )) x soi + (M 2 /(M 2 + N2 2 )y 

third process representing a new first image pri- xsoz (i) 

marily composed of the patterns of the first tis- 30 

sues of the object ^ weight factors N2 2 /(N1 2 +N2 2 ) and 

By applying the second or third method for forming N1 2 /(N1 2 +N2 2 ) in Formula (1) are determined in the 

an energy subtraction image in accordance with the manner described below. 

present invention, a new second image signal can be Specifically, if a superposition image signal S is gen- 
generated ultimately which represents a new second 35 erated from a weighted addition expressed as 
image primarily composed of the patterns of the second 

tissues of the object Specifically, the present invention s=axSOi+(\-a)xS02 p) 
still further provides a fourth method for forming an 

energy subtraction image comprising the steps of: wherein a represents an arbitrary fixed number, noise N 

i) carrying out the processes in the second or third 40 of su P*I»ation image signal S will be expressed as 
method for forming an energy subtraction image in 

accordance with the present invention, and 

ii) generating a new second image signal by carrying vi ■ (3) 

out the second process or the new second process « N (a x Nip + ((l - a) x Njp 

in which the new first image signal obtained from 45 

the third process or the new third process is taken wherein Nl represents noise contained in the first image 

as the first image signal in the second process or the signal SOI, and N2 represents noise contained in the 

new second process, the new second image signal second image signal S02. 

thus most recently generated representing a new Therefore, when the fixed number a is set such that 

second image primarily composed of the patterns 50 noise N become minimum, and the fixed number a 

of the second tissues of the object tous set is substituted into Formula (2), a superposition 

Each of the second to fourth methods for forming an image signal S can be obtained which yields a radiation 

energy subtraction image in accordance with the pres- image having good image quality. 

ent invention includes steps similar to those of the first ^ order for noise N to be minimized, the value of N 2 

method for forming an energy subtraction image in 55 m * v he minim ized Formula (3) can be transformed into 

accordance with the present invention. Therefore, as 

described above with reference to the first method for A?"(«xM)*+«i-«)xA3)? (4) 
forming an energy subtraction image in accordance 

with the present invention, the second to fourth meth- Partial differentiation of Formula (4) with respect to the 

ods for forming an energy subtraction image in accor- 60 numDer a yields 
dance with the present invention embraces various, 

substantially identical embodiments. Also, other steps, 2 2 (5) 

such as noise reducing processes, may be carried out *o ^ ~ 2(1 ~~ a ^ N1 
before or after the first to fourth methods for forming an 

energy subtraction image in accordance with the pres- 65 Noise N 2 becomes minimum a t a point, at which the 

ent invention are carried out value obtained from partial differentiation of the second 

The present invention also provides a fifth method power of noise N with respect to the fixed number a 

for forming an energy subtraction image, wherein each becomes equal to zero. Therefore, from Formula (5), a 
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value of the fixed number a is calculated which satisfies 
the condition 

2axM 2 -2<l-a)xA2 2 =0 (€) 

The following formulas obtain: 5 

a=N2 2 /(Nl 2 +N2 l ), (1 -a)=JVl 2 /(7Vl 2 +W2 2 ) (7) 

Substitution of the value of the fixed number a into 
Formula (2) yields Formula (1). Therefore, when the 10 
superposition image signal S is generated by carrying 
out the weighted addition expressed as Formula (1), 
noise of the superposition image signal S becomes mini- 
mum. 

In the fifth method for forming an energy subtraction 15 
image in accordance with the present invention, the 
superposition image signal S is generated by carrying 
out the weighted addition expressed as 



10 



5=(5JfcZ/(54H-SA2))XS01+(5Jtl/(5fcl+5Jt2))- 
XSOl 



20 



(8) 



In Formula (8), the sensitivity Ski and the sensitivity 
Sk2, which are parameters used when image processing, 
or the like, is carried out on image signals, are employed 25 
for the reasons described below. 

Specifically, in a region in which the level of the 
image signal is most predominant, the second power of 
noise has an approximately linear relationship to the 
amount of radiation, to which the stimulable phosphor 30 
sheet is exposed. Also, the sensitivity has a strong corre- 
lation to the amount of radiation, to which the stimula- 
ble phosphor sheet is exposed. Therefore, the values of 
Nl 2 and N 2 employed in the first method for forming an 
energy subtraction image in accordance with the pres- 35 
ent invention are replaced by the sensitivity Ski and the 
sensitivity Sk2, and noise contained in the superposition 
image signal S is thereby minimized. 

As indicated by line 70 illustrated in FIG. 11, the 
sensitivity, which has a linear relationship to the second ^ 
power of noise, is expressed in terms of the logarithmic 
value of the sensitivity, which is determined from the 
image signal. Therefore, the sensitivity Ski and the 
sensitivity Sk2 in Formula (8) are expressed in terms of 
the logarithmic values. Accordingly, in cases where the 45 
weighted addition is carried out by using the values of 
the sensitivity before being converted into the logarith- 
mic values, Formula (8) is transformed into the formula 



s=(iog s*r/<iog s*r+iog s*r))xs<w-t-0og 
skv/Qog skt +ios siayxsoz 



wherein each of Ski' and Sk2' represents the sensitivity 
before being converted into the logarithmic value. 

The terms "first image" and "second image" (or the 
terms "new first image" and "new second image") as 55 
used herein for the aforesaid methods for forming an 
energy subtraction image in accordance with the pres- 
ent invention mean two images, which have been ob- 
tained from energy subtraction processing and in which 
the patterns of different tissues of a single object have 60 
been emphasized or only such patterns are illustrated. 
The first image and the second image (or the new first 
image and the new second image) are not limited to 
specific images. For example, the first image and the 
second image (or the new first image and the new sec- 65 
ond image) may be a soft tissue image and a bone image. 
Alternatively, in cases where the object is a mamma of 
a human body, the first image and the second image (or 



the new first image and the new second image) may be 
an image, in which the patterns of mammary glands 
have been emphasized, and an image, in which the pat- 
tern of a malignant tumor has been emphasized. 

The present invention further provides a sixth 
method for forming an energy subtraction image, 
wherein each of the first to fourth methods for forming 
an energy subtraction image in accordance with the 
present invention is modified such that the superposi- 
tion image signal is generated by calculating a sensitiv- 
ity Ski of each of image signal components of the first 
original image signal, which represent picture elements 
in the corresponding radiation image, and a sensitivity 
Sk2 of each of image signal components of the second 
original image signal, which represent picture elements 
in the corresponding radiation image, on the basis of the 
first original image signal SOI and the second original 
image signal S02, and 
carrying out a weighted addition on the image signal 
components of the first original image signal and 
the second original image signal, which represent 
corresponding picture elements in the two kinds of 
the radiation images, the weighted addition being 
expressed as 

J=(5fc2/(5AH-5Jfc2))x5Ql+(5Jfcl/(5W+5Jk2)> 
XSOl 

The present invention still further provides a seventh 
method for forming an energy subtraction image, 
wherein each of the first to fourth methods for forming 
an energy subtraction image in accordance with the 
present invention is modified such that the superposi- 
tion image signal is generated by generating an unsharp 
signal of the first original image signal SOI and an un- 
sharp signal of the second original image signal S02, 
calculating a sensitivity Ski of each of image signal 
components of the first original image signal, 
which represent picture elements in the corre- 
sponding radiation image, and a sensitivity Sk2 of 
each of image signal components of the second 
original image signal, which represent picture ele- 
ments in the corresponding radiation image, on the 
basis of the unsharp signal of the first original 
image signal SOI and the unsharp signal of the 
second original image signal S02, and 
carrying out a weighted addition on the image signal 
components of the first original image signal and 
the second original image signal, which represent 
corresponding picture elements in the two kinds of 
the radiation images, the weighted addition being 
expressed as 

S={Sk2/{Skl+Sk2))XS<n+(SkV(jSkl+Sk2))- 
XSOl 

With the first method for forming an energy subtrac- 
tion image in accordance with the present invention, 
two kinds of radiation images of the object are formed 
with two kinds of radiation having different energy 
levels. The object is constituted of a plurality of tissues 
exhibiting different levels of radiation absorptivity with 
respect to the two kinds of radiation having different 
energy levels. The first and second original image sig- 
nals representing the two kinds of the radiation images 
are then detected. Also, the superposition image signal 
is generated by carrying out the weighted addition 
expressed as 
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„ Thereafter, a noise reducing process is carried out on 

S x$£ XS01+( * the second image signal in order to generate the second 

smoothed image signal representing a second smoothed 
wherein SOI represents the first original image signal, , ^ m which noise components of the second image 
Nl represents noise contained in the first original image have been reduced. The new first image signal is then 
signal, S02 represents the second original image signal, generated by subtracting the second smoothed image 
and N2 represents noise contained in the second original signal and the superposition image signal from each 
image signal. Thereafter, the first image signal, which other. With the two noise reducing processes, noise 
represents the first image primarily composed of pat- ir . components can be reduced in appropriate manners, 
terns of first tissues of the object, is generated by carry- Therefore, an image can be obtained which contains 
ing out the subtraction process on the first and second less noise and which has better image quality and can 
original image signals. The first smoothed image signal serve as a more effective tool in, particularly, the effi- 
is then generated by smoothing the first image signal. cient and accurate diagnosis of an illness, than the first 
The first smoothed image signal represents a first method for forming an energy subtraction image in 
smoothed image in which noise components of the first accordance with the present invention, 
image have been reduced or eliminated. Thereafter, the With the third method for forming an energy subtrac- 
first smoothed image signal is subtracted from the su- tion image in accordance with the present invention, the 
perposition image signal, and the second image signal is steps of the second method for forming an energy sub- 
thereby generated which represents the second image traction image in accordance with the present invention 
primarily composed of patterns of second tissues of the ™ are carried out repeatedly such that more noise compo- 
object Therefore, an image signal can be obtained in nents ^ ^ red uced. The respective noise reducing 
which noise components have been reduced even fur- processes can be allotted with appropriate modes of 
tner toan a superposition image signal obtained with the processing. Therefore, an image can be obtained in 
T** to f»f° r 8eneratmg a superposmon .mage signal in whkh aoi$e ^po^ ^ reduced even fur . 
the method for forming an energy subtraction image ^ T r 

SEE?! * P" 00 ** US ' appncati0 £ S £' N ^ With the fourth method for forming an energy sub- 

654,450. In this manner, a second image can be obtained . ^ . , , , , *V 

in which noise has been reduced to approximately the tractlon . m with the present inven- 

same level as that in the original images and which has twm ' aftcr ^ ste P s of 860011(1 or xhlld method for 

good image quality and can serve as an effective tool in, 30 formm S an CTCr Sy subtraction image in accordance 

particularly, the efficient and accurate diagnosis of an ^ P resent invention have been carried out, a 

illness noise reducing process is carried out on the new first 

In order for the second image having good image image signal, which has been generated by the second 

quality to be obtained, it is necessary that, in the course or tnird method for forming an energy subtraction im- 

of generating the first smoothed image signal represent- 35 age. ^ this manner, a new first smoothed image signal 

ing the first smoothed image, the signal components of & generated. Thereafter, the new first smoothed image 

the first image signal representing the patterns of the signal and the superposition image signal are subtracted 

first tissues of the object can be kept uneliminated, and from each other. Accordingly, the new second image 

only the noise components of the first image signal can can be obtained in which noise components have been 

be eliminated. However, part of the spatial frequency 40 reduced. 

components corresponding to the patterns of the first With the fifth method for forming an energy subtrac- 

tissues and part of the spatial frequency components tion image in accordance with the present invention, the 

corresponding to the noise components are identical sensitivity of the first original image signal and the sen- 

with each other. Therefore, even if a non-linear filter is sitivity of the second original image signal are detected, 

utilized which eliminates as many noise components as 45 ^ ^ and original signals m wc i g hted 

possible, the noise components and the signal compo- m accordance with the sensitivity values and added to 

nents of the first image signal representing the patterns ^ ch QthcTi md the superposition image signal is 

of the first tissues of die object cannot be completely gene rated. With the fifth method for forming 

separated from each other. an energy subtraction image in accordance with the 

Accordingly, with the second to fourth methods for 50 % ' oc • T: e ~ f . - , _ . . 

forming an Energy subtraction image in accordance V mvcntl0n ' ** m * c fmt to fourth mcthods for 
. . * " " ~T^,r . . * forming an energy subtraction image in accordance 
with the present invention, instead of aiming at com- Al _ Al f . . . 6 , , ! 
pletely eliminating the noise componenTonly with a with the present invention, an image can be obtained in 
single noise reducing process, a plurality of noise reduc- *f n ? ls ? exponents have been reduced, 
ing processes are carried out sequentially such that an 55 Wlth ™ c mth method for forming an energy subtrac- 
image can be obtained in which noise has been reduced tion una « e m accordance with the present invention, the 
and which has good image quality and can serve as an values of the sensitivities, from which the weight fac- 
efTective tool in, particularly, the efficient and accurate tors arc determined such that noise contained in the 
diagnosis of an illness. superposition image signal may be minimized, are calcu- 
Specifically, with the second method for forming an 60 lated for ^ respective image signal components of the 
energy subtraction image in accordance with the pres- first and second original image signals, which represent 
ent invention, a noise reducing process is carried out on picture elements in the corresponding radiation images, 
the first image signal in order to generate the first Therefore, noise can be minimized for each of the pic- 
smoothed image signal representing a first smoothed ture elements in the resulting image such that the image 
image in which noise components of the first image 65 quality of the image may not be adversely affected by 
have been reduced. The second image signal is then variations in the amount of radiation, which impinged 
generated by subtracting the first smoothed image sig- upon a recording medium after being irradiated to the 
nal and the superposition image signal from each other. object, for different positions on the recording medium. 



05/12/2004, EAST Version: 1.4.1 



5,402,338 

13 14 

Accordingly, the image can be obtained which has which represent radiation images of a single object, are 

good image quality. obtained, and an addition signal is obtained by adding 

With the seventh method for forming an energy sub- the image signal components of the first and second 

traction image in accordance with the present inven- image signals to each other which represent corre- 

tion, unsharp signals of the first and second original 5 sponding picture elements in the radiation images, 

unage signals are generated. The values of the sensitivi- me superposition processing method for radiation 

Ues of the respective image signal components of the images comprising generating the addition signal 

first and second original image sisals, which represent by carT ^ n out a weighted addition expressed as 
picture elements in the corresponding radiation images, 

are calculated on the basis of the unsharp signals of the io 5«(5«/(SfclH-5Jk2))x5H-(^l/(5Jtn-5A2))x52 
first and second original image signals. Therefore, the 

values of the sensitivities of the respective image signal wherein SI represents the first image signal, Ski repre- 

components of the first and second original image sig- sects a sensitivity of the first image signal, S2 represents 

nals which represent picture elements in the corre- ±c second ^ signal> Sk2 reprcscnts a sensitivity of 

sponding radiation images, can be obtained such that 15 the second image signal, and S represents the addition 

the sensitivities may not be adversely affected by minute signal 

changes and noise in the image signals. Also , appropri- By way of example, in cases where an image signal is 

ate weight factors can be determined from the calcu- a~**~*~a <v „ ui u 1. * *t_ * 

lated values of the sensitivities. SSSJJ Tf L P^hor sheet, the term 

The present invention also provides a first superposi- 20 **%J* ^ d herem „ for 

tion processing method for radiation images, wherein su P! 51 * 0Sl * 0n processmg methods and apparatuses for 

B * .• j . . . " " TT*"' , radiation images m accordance with the present inven- 

first and second image signals, each of which is made up „„„ „,„,„„ ■ - , i 

of a series of image signalcomponents and which repri- *T mMn ? ^nsitivity of the unage signal. Li cases 

sent radiation imTges of a singte object, are obtained, 811 ^ age ** cte *? *"* f^f^ ° n 

and an addition signal is obtained by adding the image « ^ ** 

^^^^Jconi^^ 25 "asusedhe^meansthe 

to each other which represent corresponding picture ' ' „ . . * y , . . „ 

elements m the radiation images, The present invention also provides an apparatus for 

the superposition processing method for radiation f^ g ,° U L * ^P^ 8 * 011 pt0 ^^ 

Jnages comprising generating ; «he addition signal 30 * ™* ** 

by carrying out a weighted addition expressed as F c .« « £ - . 

r Specifically, the present mvention also provides a 

s=<JV2V(M 2 + JV2 2 )) xsi -KJViVyvi 2 + W2 2 ))x - second superposition processing apparatus for radiation 

£2 images comprising a means for obtaining first and sec- 
ond image signals, each of which is made up of a series 

wherein SI represents the first image signal, Nl repre- 35 of image signal components and which represent radia- 

sents noise contained in the first image signal, S2 repre- tion images of a single object, and an addition means for 

sents the second image signal, N2 represents noise con- obtaining an addition signal by adding the image signal 

tained in the second image signal, and S represents the ' components of the first and second image signals to 

addition signal. each other which represent corresponding picture ele- 

The present invention further provides an apparatus ments in the radiation images, 

for carrying out the first superposition processing wherein the addition means carries out a weighted 

method for radiation images in accordance with the addition expressed as 
present invention. 

Specifically, the present invention further provides a s=(sWiSki+skX))xsi+{SkV(Ski+sk2))xS2 
first superposition processing apparatus for radiation 45 

images comprising a means for obtaining first and sec* wherein SI represents the first image signal, Ski repre- 

ond image signals, each of which is made up of a series sents a sensitivity of the first image signal, S2 represents 

of image signal components and which represent radia- the second image signal, Sk2 represents a sensitivity of 

tion images of a single object, and an addition means for the second image signal, and S represents the addition 

obtaining an addition signal by adding the image signal 50 signal. 

components of the first and second image signals to In the first superposition processing method and ap- 

each other which represent corresponding picture ele- paratus for radiation images in accordance with the 

ments in the radiation images, present invention, the addition signal S is generated 

wherein the addition means carries out a weighted from the weighted addition expressed as 
addition expressed as 55 

s=(vv2 2 /(m 2 +jvS 2 ))xsi+(AaV(Jvi2+^2 2 ))x- 

^(MV^+AZ^XSt-fCMV^+Aa^X- 52 W 
S2 

The weight factors N2V(Nl2-i-N2*) and 
wherein SI represents the first image signal, Nl repre- 60 N1 2 /(N1 2 +N2 2 ) in Formula (9) are determined in the 
sents noise contained in the first image signal, S2 repre- manner described above with reference to Formulas (2) 
sents the second image signal, N2 represents noise con- through (7). In this case, SOI and S02 in Formula (2) 
tained in the second image signal, and S represents the and in the corresponding descriptions should be re- 
addition signal. placed respectively by SI and S2, and the term "super- 
The present invention still further provides a second 65 position image signal S" in the descriptions correspond- 
superposition processing method for radiation images, ing to Formulas (2) through (7) should be replaced by 
wherein first and second image signals, each of which is the term "addition signal S." Substitution of the value of 
made up of a series of image signal components and the fixed number a, which has been obtained from For- 
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mulas (7), into Formula (2), in which SOI and S02 have 
been replaced respectively by SI and S2 ( yields For- 
mula (9). Therefore, when the addition signal S is gener- 
ated by carrying out the weighted addition expressed as 
Formula (9), noise of the addition signal S becomes 5 

minimum . 

In the second superposition processing method and 
apparatus for radiation images in accordance with the 
present invention, the addition signal S is generated by 
carrying out the weighted addition expressed as 10 

S-(SkZ/(Skl+Sk2))xSl+(Skl/(Skl+Sk2))xS2 (10) 

In Formula (10), the sensitivity Ski and the sensitivity 
Sk2, which are parameters used when image processing, i s 
or the like, is carried out on image signals, are employed 
for the reasons described below. 

Specifically, in a region in which the level of the 
image signal is most predominant, the second power of 
noise has an approximately linear relationship to the 20 
amount of radiation, to which the recording medium for 
recording the image thereon is exposed. Also, the sensi- 
tivity has a strong correlation to the amount of radia- 
tion, to which the recording medium is exposed. There- 
fore, the values of Nl 2 and N 2 employed in the first 25 
superposition processing method and apparatus for ra- 
diation images in accordance with the present invention 
are replaced by the sensitivity Ski and the sensitivity 
Sk2, and noise contained in the addition signal S is 
thereby minimized. 30 

As indicated by line 70 illustrated in FIG. 11, the 
sensitivity, which has a linear relationship to the second 
power of noise, is expressed in terms of the logarithmic 
value of the sensitivity, which is determined from the 
image signal. Therefore, the sensitivity Ski and the 35 
sensitivity Sk2 in Formula (10) are expressed in terms of 
the logarithmic values. Accordingly, in cases where the 
weighted addition is carried out by using the values of 
the sensitivity before being converted into the logarith- 
mic values, Formula (10) is transformed into the for- W 
mula 

5=(log Sfc2'/0og JAT+log 5Jk20)X«+0og 
S*lVGog S*l'+log SkiyxSZ (10-) 

43 

wherein each of Ski' and Ski' represents the sensitivity 
before being converted into the logarithmic value. 

The present invention further provides a third super- 
position processing method for radiation images, 
wherein first and second image signals, each of which is 50 
made up of a series of image signal components and 
which represent radiation images of a single object, are 
obtained, and an addition signal S is obtained by adding 
the image signal components of the first and second 
image signals to each other which represent cone- 55 
sponding picture elements in the radiation images, 
the superposition processing method for radiation 
images comprising generating the addition signal S 
by the steps of: 
i) calculating a sensitivity Ski of each of image signal 60 
components of the first image signal, which repre- 
sent picture elements in the corresponding radia- 
tion image, and a sensitivity Sk2 of each of image 
signal components of the second image signal, 
which represent picture elements in the corre- 65 
sponding radiation image, on the basis of the first 
image signal SI and the second image signal S2, 
and 
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ii) carrying out a weighted addition on the image 
signal components of the first image signal and the 
second image signal, which represent correspond- 
ing picture elements in the two kinds of the radia- 
tion images, the weighted addition being expressed 
as 

S=(Sk2/{$kl+Sk2))xSl+(Skl/(Skl+Sk2))xS2 

The present invention still further provides a third 
superposition processing apparatus for radiation images 
comprising a means for obtaining first and second image 
signals, each of which is made up of a series of image 
signal components and which represent radiation im- 
ages of a single object, and an addition means for obtain- 
ing an addition signal S by adding the image signal 
components of the first and second image signals to 
each other which represent corresponding picture ele- 
ments in the radiation images, 

wherein the addition means comprises: 

i) an operation device for calculating a sensitivity Ski 
of each of image signal components of the first 
image signal, which represent picture elements in 
the corresponding radiation image, and a sensitiv- 
ity Sk2 of each of image signal components of the 
second image signal, which represent picture ele- 
ments in the corresponding radiation image, on the 
basis of the first image signal SI and the second 
image signal S2, and 

ii) an addition device for carrying out a weighted 
addition on the image signal components of the 
first image signal and the second image signal, 
which represent corresponding picture elements in 
the two kinds of the radiation images, the weighted 
addition being expressed as 

S=@k2/(Skl+Sk2))xSi+(Skl/(Skl+Skl»XS2 

The present invention also provides a fourth superpo- 
sition processing method for radiation images, wherein 
first and second image signals, each of which is made up 
of a series of image signal components and which repre- 
sent radiation images of a single object, are obtained, 
and an addition signal S is obtained by adding the image 
signal components of the first and second image signals 
to each other which represent corresponding picture 
elements in the radiation images, 

the superposition processing method for radiation 
images comprising generating the addition signal S 
by the steps of; 

i) generating an unsharp signal of the first image sig- 
nal SI and an unsharp signal of the second image 
signal S2, 

ii) calculating a sensitivity Ski of each of image signal 
components of the first image signal, which repre- 
sent picture elements in the corresponding radia- 
tion image, and a sensitivity Sk2 of each of image 
signal components of the second image signal, 
which represent picture elements in the corre- 
sponding radiation image, on the basis of the un- 
sharp signal of the first image signal SI and the 
unsharp signal of the second image signal S2, and 

iii) carrying out a weighted addition on the image 
signal components of the first image signal and the 
second image signal, which represent correspond- 
ing picture elements in the two kinds of the radia- 
tion images, the weighted addition being expressed 
as 



05/12/2004, EAST Version: 1.4.1 



5,402,338 

17 .18 

first and second image signals. Therefore, the values of 
s=(5«/(5« + 5«))xa + (S*V(S« + s*2))xs ±e of £ s ^tive image signal compo- 
se present invention further provides a fourth su- ^l^l™*?™^ 5 ' s «»^ which 

. , y ~ " , sent picture elements in the corresponding radiation 

perposmon processing apparatus for radiation unages 5 ■ ^ ^ ^ obtained such that ^sensitivities 

comr^gamea^forobtammgfirstand second unage * ^ 

signals, each of which is made up of a series of image • . • « . - . wua,i »7 uuiac 

signal components and which represent radiation im- m * e ^P^^^ *«to» 

ages of a single object, and an addition means for obtain- * e **e™med from the calculated values of the 

ing an addition signal S by adding the image signal 10 sensmvities - 

components of the first and second image signals to BRIEF DESCRIPTION OF THE DRAWINGS 

each other which represent corresponding picture ele- . 

ments in the radiation images, mG ' } B a flow chart showing the processes, which 

wherein the addition means comprises: m carried out m a* ^age processing and displaying 

i) an operation device for generating an unsharp sig- 15 apparatus,^ 

nal of the first image signal SI and an unsharp 2 is a graph showing spatial frequency spectra 
signal of the second image signal S2, and °^ a DOne image and an image obtained by processing 
calculating a sensitivity Ski of each of image signal ^ e b° ne image signal representing the bone image, 
components of the first image signal, which repre- FIGS. 3A and 3B are graphs showing examples of 
sent picture elements in the corresponding radia- 20 probability density functions of image signal compo- 
tion image, and a sensitivity Sk2 of each of image nents, which image signal components represent a phi- 
signal components of the second image signal, rality of picture elements belonging to a predetermined 
which represent picture elements in the corre- region having a predetermined picture element in the 
spending radiation image, on the basis of the un- middle, 

sharp signal of the first image signal SI and the 25 FIG. 4 is a graph showing an example of a function, 

unsharp signal of the second image signal S2, and in which the difference between the value of an image 

ii) an addition device for carrying out a weighted signal SI and the value SI' of the image signal compo- 
addition on the image signal components of the nent representing a predetermined picture element lo- 
first image signal and the second image signal, cated in the middle of a predetermined region serves as 
which represent corresponding picture elements in a variable, 

the two kinds of the radiation images, the weighted FIG. 5 is a flow chart showing the processes, which 

addition being expressed as are carried out in a different embodiment of the method 

s~(Ski/(Ski+s!tiy>*si+(ski/(Skt 4. jjKm v « f ° r formm 6 ^ subtraction image in accordance 

S-{Sk2/(ski^sn))xsiHSki/(Ski^sk2))xS2 ^ ^ pregent mventioni 

With the superposition processing methods and appa- 35 J*»- P ttaM * ™ m f a P h ? ^ 

ratuses for ra^ation images in accordance with the f g pr0fil f * ? ^ whlch m & oym m FIG * 

present invention, the weight factors are determined 5 ' ^ng ^predetermmed ejection, 

such that noise contained in the addition signal S may be . FI . G ' 7 15 a ^f* 1 1* 0 ™** P rofilc ° f » superposi- 

minimized. The weight factors are utilized when the JL^f * Slg ^' « 

image signals for yielding the addition signal are 40 FIG. 8 is a flow chart showing the processes, which 

weighted. Therefore, noise contained in the image re- aresubstantially identical with those shown in FIG. 5, 

produced from the addition signal can be nimimized, 9 * a schematic showing an X-ray image 

and an image having good graininess and good image recording apparatus, 

quality can be obtained. FIG - *0 is a perspective view showing an X-ray 

Also, with the superposition processing methods and 45 unage read-out apparatus and an image processing and 

apparatuses for radiation images in accordance with the displaying apparatus, in which the method for forming 

present invention, the values of the sensitivities, from 311 ener £y subtraction image in accordance with the 

which the weight factors are determined such that noise present invention may be employed, 

contained in the addition signal S may be mmimized, are ^9*.** a g^P* 1 lowing the relationship between 

calculated for the respective image signal components 50 sensitivity and noise, 

of the first and second image signals, which represent FIQ. 12 is a graph showing the relationship between 

picture elements in the corresponding radiation images. sensitivity and an image signal, 

Therefore, noise can be rninimized for each of the pic- FIQ. 13 is a graph showing a stimulable phosphor 

ture elements in the resulting image such that the image sheet on an x-y coordinate system, 

quality of the image may not be adversely affected by 55 FIG. 14 is a schematic view showing how radiation 

variations in the amount of radiation, which impinged images are recorded in an embodiment of the superposi- 

upon a recording medium after being irradiated to the tion processing method for radiation images in accor- 

object, for different positions on the recording medium. dance with the present invention, 

Accordingly, the image can be obtained which has FIG. 15 is a schematic view showing how a radiation 

good image quality. 60 image is read out from a stimulable phosphor sheet, on 

With the fourth superposition processing method and which the radiation image has been stored, and 

apparatus for radiation images in accordance with the FIG. 16 is a block diagram showing an apparatus for 

present invention, unsharp signals of the first and sec- carrying out superposition processing. 

ond image signals are generated. The values of the sen- 

sitivities of the respective image signal components of 65 DESCRIPTION OF THE PREFERRED 

the first and second image signals, which represent EMBODIMENTS 

pictu/e elements in the corresponding radiation images, The present invention will hereinbelow be described 

are calculated on the basis of the unsharp signals of the in further detail with reference to the accompanying 
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drawings. First, embodiments of the method for form- 
ing an energy subtraction image in accordance with the 
present invention will be described below. In the em- 
bodiments described below, stimulable phosphor sheets 
are utilized as the recording media. 

FIG. 9 is a schematic view showing an X-ray image 
recording apparatus 1. 

With reference to FIG. 9, X-rays 3 are produced by 
an X-ray tube 2 of the X-ray image recording apparatus 
1 and irradiated to an object 4 (in this example, the chest 
of a human body). X-rays 3a, which have passed 
through the object 4, impinge upon a first stimulable 
phosphor sheet 5, and energy from the comparatively 
low energy components of the X-rays 3a is stored on 
the first stimulable phosphor sheet 5. In this manner, an 
X-ray image of the object 4 is stored on the first stimula- 
ble phosphor sheet 5. X-rays 3b, which have passed 
through the first stimulable phosphor sheet 5, then pass 
through a Miter 6 for filtering out the low energy com- 



20 
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normal to the sub-scanning direction indicated by the 
arrow Y. When the first stimulable phosphor sheet 5 is 
exposed to the laser beam 17, the exposed portion of the 
first stimulable phosphor sheet 5 emits light 22 in an 
amount proportional to the amount of energy stored 
thereon during its exposure to the X-rays. The emitted 
light 22 is guided by a light guide member 23, and pho- 
toclectrically detected by a photomultiplier 24. The 
light guide member 23 is made from a light guiding 
material, such as an acrylic plate. The light guide mem- 
ber 23 has a linear light input face 23a, positioned to 
extend along the main scanning line on the first stimula- 
ble phosphor sheet 5, and a ring-shaped light output 
face 23b, positioned so that it is in close contact with a 
light receiving face of the photomultiplier 24. The emit- 
ted light 22, which has entered the light guide member 
23 from its light input face 23a, is guided through re- 
peated total reflection inside of the light guide member 
23, emanates from the light output face 236, and is re- 



poncnts of the X-rays. X-rays 3c, which have passed 20 ceived by the photomultiplier 24. In this manner, the 



through the filter 6 and are composed of the high en- 
ergy components, impinge upon a second stimulable 
phosphor sheet 7. In this manner, an X-ray image of the 
object 4 is stored on the second stimulable phosphor 
sheet 7. During the image recording operation, marks 8, 25 
8 are placed on the object 4. The images of the marks 8, 
8 are utilized in the course of adjusting the positions of 
the two X-ray images so that the two X-ray images 
coincide with each other. 

In the X-ray image recording apparatus 1, the X-ray 
images are stored on the first stimulable phosphor sheet 
5 and the second stimulable phosphor sheet 7 with a 
single, simultaneous recording operation. Alternatively, 
the two X-ray images may be recorded one after the 
other with two independent recording operations. 

FIG. 10 is a perspective view showing an X-ray 
image read-out apparatus 10 and an image processing 
and displaying apparatus 30, in which the method for 
forming an energy subtraction image in accordance 
with the present invention may be employed. 

After the first X-ray image and the second X-ray 
image have been stored respectively on the first stimula- 
ble phosphor sheet 5 and the second stimulable phos- 
phor sheet 7 in the X-ray image recording apparatus 1 



amount of the emitted light 22, which amount repre- 
sents the first X-ray image stored on the first stimulable 
phosphor sheet 5, is converted into an electric signal by 
the photomultiplier 24. 

An analog signal S generated by the photomultiplier 
24 is logarithmically amplified by a logarithmic ampli- 
fier 25, and fed into an A/D converter 26. The A/D 
converter 26 samples the analog signal S, and the sam- 
pled signal is converted into a digital image signal SO. 
The image signal SO thus obtained represents the first 
X-ray image, which was stored on the first stimulable 
phosphor sheet 5, and will hereafter be referred to as the 
first image signal SOI (the first original image signal). 
The first image signal SOI is stored in an internal mem- 
ory of the image processing and displaying apparatus 
30. 

The image processing and displaying apparatus 30 is 
provided with a keyboard 31, from which various in- 
structions are entered, and a CRT display device 32, 
40 which displays auxiliary information for instructions 
and a visible image represented by an image signal. The 
image processing and displaying apparatus 30 is also 
provided with a floppy disk drive unit 33, which re- 
ceives and operates a floppy disk serving as an auxiliary 



30 



35 



shown in FIG. 9, the first stimulable phosphor sheet 5 45 storage medium, and a main body 34 which incorpo- 



and the second stimulable phosphor sheet 7 are placed 
one after the other at a predetermined position in the 
X-ray image read-out apparatus 10 shown in FIG. 10. 
How the first X-ray image is read out from the first 
stimulable phosphor sheet 5 will be described hereinbe- 
low. 

With reference to FIG. 10, the first stimulable phos- 
phor sheet 5 is conveyed in a sub-scanning direction 
indicated by the arrow Y by a sheet conveyance means 



rates a CPU and the internal memory. 

Thereafter, in the same manner as that described 
above, a second image signal S02 (a second original 
image signal) is obtained which represents the second 
50 X-ray image stored on the second stimulable phosphor 
sheet 7. The second image signal S02 is stored in the 
internal memory of the image processing and displaying 
apparatus 30. 

FIG. 1 is a flow chart showing the processes, which 



15, which is constituted of an endless belt, or the like, 55 are carried out in the image processing and displaying 
and which is operated by an operating means (not — 
shown). A laser beam 17, which serves as stimulating 
rays, is produced by a laser beam source 16. The laser 
beam 17 is reflected and deflected by a rotating polygon 
mirror 19, which is being quickly rotated by a motor 18 60 
in the direction indicated by the arrow Z. The laser 
beam 17 then passes through a converging lens 20, 
which may be constituted of an f$ lens, or the like. The 
direction of the optical path of the laser beam 17 is then 
changed by a mirror 21, and the laser beam 17 is caused 65 
to impinge upon the first stimulable phosphor sheet 5 
and scan it in a main scanning direction indicated by the 
arrow X. The main scanning direction is approximately 



apparatus 30. The processes are carried out on the first 
image signal SOI representing the first X-ray image and 
the second image signal S02 representing the second 
X-ray image, which signals are stored in the internal 
memory of the image processing and displaying appara- 
tus 30. 

The first image signal SOI and the second image 
signal S02, which are stored in the internal memory of 
the image processing and displaying apparatus 30, rep- 
resent a first X-ray image 41 and a second X-ray image 
42 shown in FIG. L The first X-ray image 41 has been 
recorded with the comparatively low energy compo- 
nents of the X-rays. The second X-ray image 42 has 
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been recorded with the comparatively high energy calculated. The sensitivity corresponds to the photoe- 
componentsofthe X-rays. Both of the first X-ray image lectric conversion factor, which represents to what 
41 and the second X-ray image 42 are original images image signal level a predetermined amount of emitted 
composed of patterns of soft tissues and bones. The light is to be converted. The sensitivity serves as a pa- 
levels of image density of the soft tissue patterns and the 5 rameter for determining the read-out conditions, under 
bone Patterns are different between the first X-ray wbich ^ ^ figDa k detected( Md ^ ^ 
image 41 and the second X-ray image 42. cessing conditions, under which an image signal is to be 

*m .Tf hT SOI and the second image processed . Operations for calculating the vllues of the 
signal S02 are read from the internal memory of the rea d-out conditions and/or the image processing condi- 
m^ge prooessmg and displaying apparatus 30 shown in m tions „ cuM m b a « 

FIG. 10. Position adjustment processing is then earned „„«u, • . : A * p * . 

out on the first image signal SOI and the second image ^Z^Tf* * 8 f ' A ^ n ^ ° f ™S e 
signal S02 such that ^positions of the first X-riy ^ ?T< n ^number of radiation images 

image 41 represented by trTfirst image signal SOI and ™ Processed. The algonthms which calcu- 

the second X-ray image 42 represented byfce second 15 !* te * e and/or the image process- 

image signal S02 may coincide with each other. For "g ?° n ff ons " on i <he **** of results 

this purpose, a method disclosed in, for example, U.S. o^ed from this processing. As examples of the algo- 
Pat No. 4,710,875 may be employed. With the position nthms whlch have heretofore ^en employed, methods 
adjustment processing, one of the two X-ray images is ^ oyni wherein a probability density function of an 

linearly moved or rotated with respect to the other 20 m * g * . SlgnaI * crcatcd and ^yzed, and the read-out 
X-ray image until the images of the marks 8, 8 in one conditions and/or the image processing conditions are 
X-ray image, which marks are shown in FIG. 9, overlap determined on the basis of the results of the analysis of 
the images of the marks 8, 8 in the other X-ray image. probability density function. Such examples of the 

Thereafter, a subtraction process is carried out on the algorithms are disclosed in, for example, U.S. Pat. No. 
first image signal SOI and the second image signal S02. 25 4,638,162 and Japanese Unexamined Patent Publication 

Specifically, X-ray absorption coefficients ji are clas- No * 61(1986>280163. Therefore, how the sensitivity 
sified into the following: can be calculated will not herein be described in detail. 

fiLT: Absorption coefficient of soft tissues with re- As indicated by line 70 shown in FIG. 11, the sensi- 
spect to the low energy components of X-rays. tivity Ski and the sensitivity Sk2 t which are calculated 

fiHT: Absorption coefficient of soft tissues with re- 3 q in this embodiment, are expressed in terms of the loga- 
spect to the high energy components of X-rays. rithmic values of the sensitivity Ski' and the sensitivity 

fiLB: Absorption coefficient of bones with respect to Sk2' detected from the image signals. Specifically, the 
the low energy components of X-rays. following formulas obtain. 

/xHB: Absorption coefficient of bones with respect to 
the high energy components of X-rays. 35 5*1 « log SfcT, 

The first image signal SOI and the second image 
signal S02 are weighted, and the image signal compo- SA2=log Sk2' 

nents of the weighted image signals are subtracted from 

each other which represent the image information After sensitivity Ski and the sensitivity Sk2 are 
stored at corresponding picture elements in the two aq calculated, the first image signal SOI and the second 
X-ray images. Thus a bone image signal SI is obtained, vaa&ge signal S02 are weighted with the sensitivity Ski 
which can be expressed as 311(1 ^ sensitivity Sk2. The image signal components of 

the weighted first and second image signals are then 
i^lj /jjv added to each other which represent corresponding 

Sl => SOl - pjyy- SOl + C 45 picture elements in the two X-ray images. The weighted 

addition is represented by the formula 
where C denotes a bias component. The bone image _ 

signal SI represents a bone image 43 shown in FIG 1, so=wW+smxsoi«skv { Ski + sk2)). 
which image is composed of the bone patterns. 1 

• ^ J£? ^ "?J** f eCOnd TP 50 In this manner, a superposition image signal SO, which 

S02 may be weighted m a duTerent way, and the represents a superposition image 44, is obtained from 

^vLTs? ZTr™ °w h Wei ^ "ff" the two image signals SOI and S02. Formula <8") cor- 
may be subtracted fromeach other which represent the ds to S For ^ a m which Nl2 is Ie ± c ^ b 

Zk ^?^ 0n v at COrr f I T dmg piCtUle de fi thesensitivity Ski, and N» is replaced by the sensitivity 
ments m the two X-ray images. In this manner, a soft 55 C i,<> m. / • * • \. • r 

tissue image signal S2 can be obtained, which is «- f^^ 0 "-' nolse H cont ^f d m *e superposition 
pressed as image 44 is minimized, and the superposition image 



signal SO is advantageous for the subsequent processes. 
In cases where the weighted addition is carried out 
sz = soi - sol + c (12) by using the values of the sensitivities Ski' and Sk2' 

^ 60 before being converted into the sensitivities Ski and 



where C denotes a bias component The soft tissue Sk2 logarithmic values, the superposition 

image signal S2 represents a softtissue image composed mage ^ S ° may * generated Wth to fonnu3a 

of the soft tissue patterns. However, in this embodi- ^-(logStt70ogS*r+log Sfa*))xSOi+Qo& 

ment, the operations for generating the soft tissue image 65 Skt'/Qog Ski' +log Skz 1 )) xS02 (n> 
signal S2 need not be carried out 

Also, a sensitivity Ski of the first image signal SOl In cases where graduations are used for the sake of 

and a sensitivity Sk2 of the second image signal S02 are convenience for the sensitivity Sk shown in FIG. 11 
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(e.g., in cases where graduations on the horizontal axis filter adaptive to a probability density function is uti- 

Sk" shown in FIG. 11 are used), the relationship be- lized. With the smoothing process, noise can be elimi- 

tween Sk" and Sk' is expressed as nated such that edges (Le. step-like changes in density, 

lxkf , which define boundaries among patterns of a plurality 

£f=4xifr 5 of different tissues of an object), which it is necessary to 

r , iV . reproduce, may be kept sharp and no artifact may occur 

Therefore such ^ Ae su^rpoation image sigiul jn ^ ^ ^ noise ^ ^ eliminated 

SO may be generated with the formula quickly ^ 

so=(4- 1 SJkry((4- ) Sjfci")+(4-5Jk2 w ))- Specifically, first, each of the picture elements in the 

xsoi+n-skvya*-, s*i")+(4-SA2"))xS02 (14) 10 t>°ne image 43 is taken as a predetermined picture ele- 
ment, and the probability density function of the image 

Thereafter, the bone image signal SI is processed signal components of the bone image signal SI is gener- 

such that noise components contained in the bone image at ed, which represent a plurality of the picture elements 

43 may be eliminated. belonging to a predetermined region having the prede- 

FIG. 2 is a graph showing spectra of a bone image 15 termined picture element in the middle, 

and an image, which is obtained by processing the bone FIGS. 3A and 3B are graphs showing examples of 

image signal representing the bone image, with respect probability density functions of image signal compo- 

to the spatial frequency, f. nents of the image signal SI, which image signal compo- 

In FIG. 2, curve 51 indicates the spectrum of the nents represent a plurality of picture elements belonging 
bone image 43, and curve 53 indicates the spectrum of 20 to a predetermined region having a predetermined pic- 
noise components included in the bone image 43. ture element in the middle. The image signal component 

First, a smoothing process is carried out on the bone representing the predetermined picture element has a 

image signal SL As the smoothing process, one of vari- value SI'. FIG. 4 is a graph showing an example of a 

ous processes may be employed. For example, a simple function, in which the difference between the value of 

averaging process may be employed wherein the mean 25 the image signal SI and the value SI' of the image signal 

value of the values of the image signal components of an component representing the predetermined picture ele- 

image signal, which represent the picture elements be- ment located in the middle of the predetermined region 

longing to a predetermined region having a predeter- serves as a variable. 

mined picture element in the middle, is calculated and The probability density functions shown in FIGS. 3A 
employed as the value of the image signal component 30 and 3B m denoted by h(Sl). Also, a function, the value 
representing the predetermined picture element. Alter- of which decreases monotonously as the absolute value 
natively, a median filter process may be employed |S1-S1'| increases, e.g. the function shown in FIG. 4, is 
wherein the median value of the values of the image denoted by f(Sl-Sl'). The values of a function g(Sl) 
signal components of an image signal, which represent representing how frequently the values of image signal 
the picture elements belonging to a predetermined re- 35 components of an image signal occur, which image 
gion having a predetermined picture element in the signal has been processed, are calculated with the for- 
middle, is calculated and employed as the value of the mu ] a 
image signal component representing the predeter- 
mined picture element As another alternative, an edge g(Si) = h(si) xfisi-si') (15) 
keeping filter (V-filter) process may be employed 40 

wherein a predetermined region having a predeter- In cases where the function h(Sl) includes a plurality of 

mined picture element in the middle is divided into a projecting parts shown in FIG. 3A, the function g(Sl) 

plurality of small regions, and the variance of the values has the effects of extracting only of the projecting part, 

of the image signal components corresponding to each to which the image signal component having the value 

small region is calculated. A small region associated 45 of SI' and representing the predetermined picture ele- 

with the smallest variance is then found, and the mean ment belongs. 

value of the values of the image signal components After the values of the function g(Sl) have been cal- 

corresponding to the small region associated with the culated with Formula (15), the values of the image 

smallest variance is employed as the value of the image signal components of the image signal SI, which image 

signal component representing the predetermined pic- 50 signal components represent the picture elements be- 

ture element As a further alternative, a process may be longing to the predetermined region, are weighted with 

employed wherein Fourier transformation is carried out the values of the function g(Sl). A calculation is then 

on an image signal, the signal obtained from the Fourier made to find a mean-level value 'Sim of the weighted 

transformation is subjected to an operation for remov- values of the image signal components of the image 

ing high spatial frequency components corresponding 55 signal SI. Specifically, by way of example, the moment 

to noise components, and thereafter inverse Fourier of first order of the function g(Sl) is calculated with the 

transformation is carried out. formula 

However, the simple averaging process has the draw- 
backs in that edges in the image become unsharp. The Sim=(Si)xsidSi/fsidSi (16) 
median filter process has the drawbacks in that, because 60 

picture elements are interchanged, contour line-like The picture elements in the bone image 43 are sequen- 

artifacts often occur. The edge keeping filter process tially token as the predetermined picture element, and 

has the drawbacks in that honeycomb-like artifacts the processes with Formulas (15) and (16) are carried 

often occur. The Fourier transformation process has the out f° r *H of the picture elements in the bone image 43. 

drawbacks in that a long time is taken for operations to 65 ^ n this manner, a smoothed image signal Sim is gener- 

be carried out. ated. (As an aid in facilitating the explanation, the same 

Therefore, in this embodiment, as will be described reference numeral is utilized to indicate both the value 

below, a smoothing process is carried out in which a °f the image signal component representing each pic- 
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ture element and the image signal representing the distributed, which image signal components represent 
whole image.) As indicated by curve 52 in FIG. 2, the the picture elements located along a predetermined 
smoothed image signal Sim is generated by primarily direction (x direction) in the first X-ray image 41. FIG. 
eliminating the high spatial frequency components from 6B shows how the values of the image signal compo- 
the bone image signal SI. As shown in FIG. 3A, as for 5 nents of the second image signal S02 representing the 
a picture ^element located I in the vicmity of an edge, the seC0 nd X -ray image 42 are distributed, which image 
smoothed image signal Sim has the mean-level value of signal components represent the picture elements lo- 
*e values belonging only to the projecting part, to ^ ^ ^ p^^ed ^ecHon (x direction) in 
which said picture element belongs. Therefore, edges in ^ second x _ ra £ ^ 42 The leyels of ^ ^ ^ 
thebone unage 43 can be kept sharp. 10 sigDal sol and die second image signal S02 are tffo! 

Thereafter, the superposition image signal SO, which ^jTZL^ m . „ B °~ e TT 1 r\ . 

is expressed as Formula (8), Formula oSTor Formula ^STnT?^ h ' ° ?SrS^ ^ 

(14) and represents the superposition image 44, and the S ?l ™ 6 .""J* ^ sl ^ S ° 2 is com- 

smoothed image signal Sim are weighted. The image ^ unage Slgn ^ ™ m P° nen .^ which represent 

signal components of the weighted smoothed image i 5 ^V 0 * Patterns ^rrespondmg to the hatched 
signal Sim are subtracted from the image signal compo- re ^ on m m ?' 6A or FIG ' m ) have approximately 
nents of the weighted superposition image signal SO, umform v f ues ' ^ components, which 

which image signal components represent the image re P resent toe bone patterns and have values changing 
information stored at corresponding picture elements in step-wise, and the random noise components. These 
the two X-ray images. Thus a soft tissue image signal 20 ***** of ima S e si S nal components are superposed 
S2' is obtained, which can be expressed as one u P on another. 

By carrying out the weighted subtraction process 
with Formula (12) on the first image signal SOI repre- 
IlLB \ senting the first X-ray image 41 (the original image) and 

( _ J + fiHB ) 25 the second image signal S02 representing the second 

52 =5 ° T^Zb uLT \ Slm * C ?" ray " nage 42 ^ e crista 1 image), the soft tissue 

l - [~{S5 f3rr~ J image signal S2 representing the soft tissue image 47 is 

generated. Also, by carrying out the addition process 

„ . ^ ■* . . „ & , with Formula (8")> Formula (13), or Formula (14) on 
where C denotes a bias component The soft tissue *u A : n em , V _. . v . , 

image signal S2< represents a proved soft tissue image 30 f^? h ™* e "J* SO } * e ™f* ma * e ^ 
46 shown in FIG. L The processed soft tissue image 46 f° 2 ' the position image signal SO representing 
has approximately the same image information £ the . . 

soft tissue image expressed as Formula (12) and includes n °' 7 • how * how v .^ ues . of me """^ ^al 
less noise components than the soft tissue image ex- ? m ? u on ^ te T °! ^e superposiUon image signal SO are 
pressed as Formula (12). 35 attributed. the first image signal SOI shown in 

The soft tissue image signal S2' which has been gener- ™' ^ second S02 shown in 

ated with Formula (17), is fed into the CRT display mG ' ™» ^ su Perposition image signal SO is corn- 
device 32 of the image processing and displaying appa- P 08 ^ of ^ "S 1131 components, which represent 
ratus 30. A visible image is reproduced from the soft ^e soft tissue patterns (corresponding to the hatched 
tissue image signal S2' and displayed on the CRT dis- 40 re gion ^ 7) and have approximately uniform val- 
play device 32. ues » ^ e image signal components, which represent the 

In the embodiment described above, the soft tissue ^° ne Patterns and have values changing step-wise, and 
image signal S2 is generated by smoothing the bone ^ e ran ? om noise components. These three types of 
image signal SI and subtracting the smoothed signal unage signal components are superposed one upon an- 
from the original image signal representing the original 45 other. However, the superposition image signal SO 
image. In cases where a bone image is to be reproduced, includes less noise components than the first image 
the soft tissue image signal S2 is generated with For- signal SOI shown in FIG. 6A and the second image 
muia (12) and then smoothed. The smoothed signal is si S n al S02 shown in FIG. 6B. 
then subtracted from the original image signal repre- FIG. €D shows how the values of the image signal 
senting the original image. In this manner, a bone image 50 components of the soft tissue image signal S2, which has 
in which noise components have been reduced can be been generated with Formula (12), are distributed. The 
obtained. soft tissue image signal S2 is primarily composed of the 

A different embodiment of the method for forming an image signal components, which represent the soft tis- 
cnergy subtraction image in accordance with the pres- sue patterns and have approximately uniform values, 
ent invention will be described hereinbelow. 55 However, the soft tissue image signal S2 includes more 

FIG. 5 is a flow chart showing the processes in a random noise components than the first image signal 
different embodiment of the method for forming an SOI shown in FIG. 6A and the second image signal 
energy subtraction image in accordance with the pres- S02 shown in FIG. 6B. 

ent invention. FIGS. 6A, 6B, 6D through 6M and FIG. FIG. 6E shows how the values of the image signal 
7 are graphs showing the profiles of the images, which 60 components of the bone image signal SI, which may be 
are shown in FIG. 5, along a predetermined direction. generated with Formula (11), are distributed. (In this 

In FIG. 5, similar elements are numbered with the embodiment, the bone image signal SI need not be 
same reference numerals with respect to FIG. 1 generated.) The bone image signal SI is primarily com- 

FIGS. 6A and 6B show the profiles of the first X-ray posed of the image signal components, which represent 
image 41 and the second X-ray image 42, which are 65 the bone patterns and have values changing step-wise, 
original images. Specifically, FIG. 6A shows how the However, like the soft tissue image signal S2 shown in 
values of the image signal components of the first image FIG. 6D, the bone image signal SI includes more ran- 
signal SOI representing the first X-ray image 41 are dom noise components than the first image signal SOI 
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shown in FIG. 6A and the second image signal S02 rising part of the smoothed bone image signal Sim' 
shown in FIG. 6B. shown in FIG. 61 is unsharp, the information represent- 
As illustrated in FIG. 5, a smoothing process 51 is ing the corresponding part of the bone image is included 
carried out on the soft tissue image signal S2, which as noise in the soft tissue image signal S2\ However, the 
represents the soft tissue image 47 and is distributed as 5 level of random noise and the level of the information, 
shown in FIG. 6D. From the smoothing process 51, a which represents the bone image and constitutes noise, 
smoothed soft tissue image signal S2m is obtained, are very low. Therefore, a series of the processes may 
which represents a smoothed soft tissue image 61 and is be finished in this step. The soft tissue image signal S2' 
distributed in the pattern shown in FIG. 6F. With the may be fed into the CRT display device 32 of the image 
smoothing process 51, the spatial frequency compo- 10 processing and displaying apparatus 30 shown in FIG. 
nents higher than a frequency of, for example, 1.0 cy- 10 , and a visible image may be reproduced from the soft 
cle/mm are eliminated from the soft tissue image signal S2 ' and displayed on the CRT dis- 
S2 representing the soft tissue image 47. p j ay device 32 

Thereafter, the superposition image signal SO and the However, in this embodiment, the same processes as 
smoothed soft tissue image signal S2m are weighted, 15 those described above are repeated even further such 
and *e weighted smoothed soft tissue image signal S2m ^ m - ^ better image quality may be ob- 
is subtracted from the weighted superposition image tained 

3°' v mma : a ^ ^ Sff iS Ate «* ***** *»g e S2' representing the 

f^p^^r^ a W S' ,f ,n "* tisSUe ™& « has *~ generated; a smoothing 

rated in FIG. <5G, the bone image signal SI' includes 20 s 53 ^ ^ Qn t[ ; * 

less random noise components than the bone image ^ n 5. j & 

1 ci u ■ ct/*^ u lt « • S2 . From the smoothing process 53, a smoothed soft 

signal SI shown in FIG. 6E. However, the bone image ; „ ■ , c - , r *\ . . . ' . . 

• 1 o^/ 1 • 1 j 1 * tissue image signal S2m is obtained, which represents a 

signal SI shghtly includes the high spatial frequency „ , u A * ^_ , . ll 

components of L soft tissue imTge47 dueto the ^mr^f ^ " d "»*7 rted m 

smoothing process carried out on the soft tissue image 25 PfT 1 » 6J. With the smoothing process 

47 53, the spatial frequency components higher than a 

A smoothing process 52 is then carried out on the of, for example, 1.5 cycle/mm are eliminated 
bone image signal SI', which has been generated in the fro £ thc i oft S2 . _ 
manner described above. With the smoothing process Thereafter the superposition image signal SO and the 
52, patterns having low contrast and falling within the 30 smoo ^ ed soft tissue image signal S2m are weighted, 
spatial frequency region higher than, for example, 0.5 * e weighted smoothed soft tissue image signal 
cycle/mm are eliminated from the bone image 62 (i.e. ? 2m K . subtracted from Ac weighted superposition 
small changes in the bone image signal SI' are elimi- sl S^ S0 * ^ naanner, a bone image signal 
nated). For this purpose, by way of example, the bone S1 " ™ obtained, which represents a bone image 66. As 
image signal SI' may be processed with a filter de- 35 illustrated m FlG - ^ne image signal SI" in- 
scribed below. Specifically, a window having an area eludes less random noise components and less informa- 
corresponding to 0.5 cycle/mm is determined for a tioil » which represents the soft tissue image and consti- 
predetermined picture element PO. From the image mtes noise » than the bone image signal SI' shown in 
signal components of the bone image signal SI' repre- 60 ' 111 045418 where a bone image is to be repre- 
senting the picture elements belonging to the window, 40 duce 4 a visible image may be reproduced from the 
the image signal components are then found the values 00116 image signal SI" and reproduced on the CRT 
of which fall within the range of: display device 32. 

the value of an image signal component Slo' repre- In this embodiment, a smoothing process 54 is then 

senting the predetennined picture element P0±a carried out on the bone image signal SI", which has 

predetermined value. 45 be* 11 generated in the manner described above. From 

Thereafter, the mean value of the image signal compo- the smoothing process 54, a smoothed bone image sig- 

nents, which have thus been found, is calculated and nal Sim" is obtained, which represents a smoothed bone 

employed as the value of a new image signal component image 67 and is distributed in the pattern shown in FIG. 

Slo' representing the predetermined picture element ^M- With the smoothing process 54, patterns having 

P0. With the smoothing process 52, a smoothed bone 50 low Contrast and falling within the spatial frequency 

image signal Sim' is obtained which represents a region higher than, for example, 1.0 cycle/mm are elim- 

smoothed bone image 63. As illustrated in FIG. 61, in inated from the bone image 66. 

the smoothed bone image signal Sim', the noise compo- Thereafter, the superposition image signal SO and the 

nents have been reduced. Also, the high spatial fre- smoothed bone image signal Sim" are weighted, and 

quency components of the soft tissue image 47 due to 53 weighted smoothed bone'image signal Sim" is sub- 

the smoothing process carried out on the soft tissue tracted from the weighted superposition image signal 

image 47 have been reduced. However, the rising part SO. In this manner, a soft tissue image signal S2" is 

of the smoothed bone image signal Sim' becomes un- obtained. As illustrated in FIG. 6L, the soft tissue image 

sharp. signal S2" includes less random noise components and 

Thereafter, the superposition image signal SO and the 60 less information, which represents the bone image and 

smoothed bone image signal Sim' are weighted, and constitutes noise, than the soft tissue image signal S2' 

weighted smoothed bone image signal Sim' is sub- shown in FIG. 6R 

tracted from the weighted superposition image signal In the manner described above, the smoothing pro- 

SO. In this manner, a soft tissue image signal S2' is cesses and the weighted subtraction processes with 

obtained which represents a soft tissue image 64. As 65 respect to the superposition image signal SO are carried 

illustrated in FIG. 6H, the soft tissue image signal S2' out repeatedly such that the bone images and soft tissue 

includes less noise components than the soft tissue images, in which noise has been reduced sequentially, 

image signal S2 shown in FIG. 6D. Also, because the may be obtained alternately. 
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FIG. 8 is a flow chart showing the processes in a still 0*) 
further embodiment, which processes are substantially 

identical with those shown in FIG. 5. In FIG. 8, similar s , (( y^LB \ f J^_\ \ 

elements are numbered with the same reference numer- \V + ^ J sl + \ l + yJiT ) SN / 

als with respect to FIG. 5. 5 
In the embodiment of FIG. 8, the series of the pro- / - r \ 

cesses for generating the bone image 62 in the embodi- 2 • ( Wr~) 

ment of FIG. 5 (which processes correspond to the 

processes described above with reference to FIG. 1, tl^„ . 

except that the bone image and the soft tissue image are 10 ' ^ O0inpOIientS Can be reduced evcn 

interchanged with each other) are replaced by the pro- Tk™„4v„ ^ A , 
cesses described below. P * P ^ere^r the same processes as those shown m 

e „ . . . rlG. 5 are earned out. 

Specifically, the bone : image signal SI representing ^ ^ embodiment of FIG. 8, only the initial pro 
the bone unage 43 and the soft tissue image signal S2 cesses in the embodiment of FIG. 5 are replaced by the 
representing a soft tissue image 47 are generated by processes described above. Such replacement may be 
canying out calculations with Formulas (11) and (12) carried out at an arbitrary stage of the processes carried 
from the first image signal SOI representing the first out repeatedly. Such embodiments are substantially 
X-ray image 41 and the second image signal S02 repre- identical with the embodiment of FIG. 5. The method 
senting the second X-ray image 42. 20 for forming an energy subtraction image in accordance 

Thereafter, in the same manner as that in the embodi- with the present invention embraces various such, sub- 
ment described above, the smoothed image signal S2m stantially identical embodiments wherein at least one 
representing the smoothed soft tissue image, in which stage of the processes is modified, 
the noise components included in the soft tissue image In the aforesaid embodiments of the method for form- 
47 have been reduced, is generated by processing the 25 ing an energy subtraction image in accordance with the 
soft tissue image signal S2 in accordance with Formulas present invention, a soft tissue image or a bone image is 
(15) and (16). The image signal components of the formed from X-ray images of the chest of a human 
smoothed image signal S2m are then subtracted from body. However, the method for forming an energy 
the image signal components of the soft tissue image subtraction image in accordance with the present inven- 
signal S2, which image signal components represent the 30 tion is not limited to the formation of the soft tissue 
image information stored at corresponding picture ele- image or the bone image, but is applicable widely when 
ments in the two X-ray images. In this manner, a noise either one or both of two images are to be obtained, in 
signal SN representing a noise image 60, which is com- which the patterns of two different tissues of a single 
posed of only the noise components, is obtained. The object have been emphasized or only such patterns are 
noise signal SN can be expressed as 35 illustrated. For example, two such images may be an 

image, in which the patterns of mammary glands have 
SN=S2-sim (18) been emphasized, and an image, in which the pattern of 

a malignant tumor has been emphasized. 
As indicated by curve 53 in FIG. 2, the noise signal SN Also, in the aforesaid embodiments of the method for 
is composed of only the noise components included in 40 forming an energy subtraction image in accordance 
the soft tissue image 47. In the smoothed image signal with the present invention, stimulable phosphor sheets 
S2m, even if the information representing the edges in are used - However, the method for forming an energy 
the soft tissue image 47 has a level of spatial frequency subtraction image in accordance with the present inven- 
as high as that of the noise components, the information tion ^ so applicable when other recording media, such 
representing the edges will not be lost. Therefore, by 45 88 x * ra y fifin (ordinarily combined with intensifying 
carrying out the calculations with Formula (18) to find screens), are used. 

the difference between the soft tissue image signal S2 Additionally, in the aforesaid embodiments of the 
and the smoothed image signal S2m, the noise signal SN me *hod for forming an energy subtraction image in 
can be obtained in which the information representing accordance with the present invention, the weight fac- 
the edges has been completely canceled. Accordingly, 50 ^ *? J* ^ durm S the weighted addition are deter- 
the noise signal SN more accurately represents only the JSff 1 * om s ^ of the first image signal 

noise components of the soft tissue image 47 than when l°l ™ •f"*™* S ^ of the second image signal 
a smoothing process was carried out such that the infor- f O^Alternatively, noise Nl and noise N2 contained in 
mation representing the edges may be lost „ on f»* ™*UJ*t- ? 01 ^J 02 may * 

Thereafter, the noise signal SN and the bone image 55 k ^ ' the weighted addmon may be carried out 
signal SI representing the tale image 43 are wdJK pomuk ST * " ^ 

I^fi? 1 ! "TT ° f ^ v Cighted ^ Furthcr > * * e embodiments of the method 
signals are added to each other, which image signal for forming m subtraction image in accordance 
components represent the miage information stored at m ^ ^ prcsent £ vention> when £ supe ^£* 
correspondmg picture elements in the two images. In image signal SO is to be obtained, the weight factors are 
this manner, a bone image signal SI' is obtained, which determined from the sermtivities Ski and Sk2 of the 
represents a processed bone image 62. The processed two original image signals SOI and S02, which repre- 
bone image 62 has approximately the same image infor- sent image information stored in the whole areas of the 
mation as the bone image 43 and includes less noise 65 corresponding radiation images. Alternatively the sen- 
components than the bone image 43. In this embodi- sitivities Ski and Sk2 of the two original image signals 
ment, the weighted additions are carried out with the SOI and S02 may be detected for the respective image 
formula signal components of the two original image signals 
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SOI and S02 which represent image information stored sharp signals of the first and second original image sig- 

at respective picture elements in the corresponding nals. In such cases, the values of the sensitivities of the 

radiation images. How the sensitivity Ski of each of the respective image signal components of the first and 

image signal components of the original image signal second original image signals, which represent picture 

SOI, which represent image information stored at re- 5 elements in the corresponding radiation images, can be 

spective picture elements in the corresponding radiation obtained such that the sensitivities may not be adversely 

image is ^calculated will be desenbed herein^low affected by minute changes and noise in the image sig- 

First, the senauvny ^Skl and a latitude Gpl ^of the ^ ^ ^ si ^ ^ ^ obtained by J^f 

first image signal SOI are calculated from the first o % -j\„ ~~ A * 6 

imagelgnal SoTwkh the method disclosed in U.S. Pat. 10 °*° nC ° f ? m0m Voces***, such asan 

No 4,638,162. The latitude corresponds to the ratio of ™ shar P "»* * median ffltcr P rocess > ^ 

the largest amount of emitted light, which is capable of fU * r P™? 88 ' f" 1 » v prOCCSS 

being accurately converted into an image signal during e Em ^ dHnents of the superposition processing method 

the image readout operation, to the smallest amount of f or ra ^ ti0 ^ m accordance with the present 

emitted light, which is capable of being accuratelycon- 15 mvcntlon ^ bc dcscnbed neronbelow. In the embodi- 

verted into an image signal during the image read-out ments described below, stimulable phosphor sheets are 

operation. For example, in cases where the sensitivity }3Scd 35 recording media. 

Ski is calculated as being 2 and the latitude Gpl is F1G - 14 shows how radiation 102, which has passed 

calculated as being 2, the sensitivity Ski and the image through a single object 101, is irradiated to two stimula- 
signal SI have the relationship shown in FIG. 12. In 2 o ble pkosP 001 " sheets 104A and 104B. 

cases where the picture element values are 10 bits, Specifically, the first stimulable phosphor sheet 104A 

Smin=0, and Smax =1,023. When the x-y coordinates and the second stimulable phosphor sheet 1MB are 

are set on the first stimulable phosphor sheet 5 as shown superposed one upon the other, and a radiation source 

in FIG. 13, the sensitivity Skl(x,y) at an arbitrary point 103 is activated to produce the radiation 102. The radia- 

having the coordinates (x,y) on the first stimulable 2 s tion 102, which has been produced by the radiation 

phosphor sheet 5 can be calculated with the formula source 103, passes through the object 101, The radiation 

102, which has passed through the object 101, impinges 

So\(x,y) - Smax + Smin u P° n *** e ^ rst stimulable phosphor sheet 104A and the 

Sk\(x,y) Smax _ Sf ^ n x + Ski second stimulable phosphor sheet 104B. In this manner, 

p 30 radiation linages of the object 101 are stored on the first 

wherein S01(x,y) represents the value of the image ^uUble phosphor sheet 104A and the second stimula- 

signal component of the first image signal SOI, which Wepnosphor sheet 1WB. 

image signal component corresponds to an arbitrary _ ^e*™** * « radiation images are read out from the 
point having the coordinates (x,y) on the first stimulable nrststimuiable phosphor sheet 104A and the second 
phosphor sheet 5. 35 stimulable phosphor sheet 104B by using an image read- 
In the manner described above, the sensitivity out means shown in FIG. 15, and image signals repre- 
Skl(x,y) corresponding to each of the picture elements senting the radiation images are thereby obtained. Spe- 
in the radiation image, which was stored on the first cifically, the first stimulable phosphor sheet 104A is 
stimulable phosphor sheet 5, is calculated with Formula moved by a sub-scanning means 109 in the sub-scanning 
(20). Thereafter, the sensitivity Sk2(x,y) corresponding 40 direction indicated by the arrow Y. The sub-scanning 
to each of the picture elements in the radiation image, means 109 may be constituted of an endless belt, or the 
which was stored on the second stimulable phosphor hxe. At the same time, a laser beam 111, which serves as 
sheet 7, is calculated from the second image signal S02. stimulating rays, is produced by a laser beam source 
The superposition image signal SO is then generated by 110. The laser beam 111 is deflected by a inning mir- 
carrying out the calculations expressed as 45 ror 112 and caused to scan the stimulable phosphor 

sheet 104A in the main scanning directions indicated by 

SC 5&? SSJ^^^S^tfSS^ ~ ™ «« the double-headed arrow X. When the stimulable phos- 

xsoi+sn^yyiski^s^xsoi (21) phor shect 1D4A is exposed t0 ^ lascr ^ ^ it 

In the manner described above, the values of the 50 ™* ^ * ™ !? ? t0 ^° n t0 *f ^ount of energy 
sensitivities, from which the weight factors are deter- !? red *™f lts ^ t0 102 ' 
mined, are calculated for the respective image signal ™? fitted hght "3 enters a light guide member 114, 
components of the first and second original image sig- wluch 18 made ^ m a transparent acrylic plate, from its 
nals, which represent picture elements in the corre- one face ' ^ ****** light 113 is guided through 
sponding radiation images. Therefore, noise can be 55 rc P«atcd total reflection inside of the light guide mem- 
minimized for each of the picture elements in the result- 114 811(1 detected by a photomultiplier 115. The 
ing image such that the image quality of the image may photomultiplier 115 generates an image signal S A cor- 
not be adversely affected by variations in the amount of responding to the amount of the emitted Ught 113, i.e. 
radiation, which impinged upon the recording medium, representing the radiation image stored on the stimula- 
such as the stimulable phosphor sheet, after being irradi- 60 b* c phosphor sheet 104A. 

ated to the object, for different positions on the record- The image signal SA is logarithmically amplified by a 

ing medium. logarithmic amplifier 116 and is then converted by an 

Also, unsharp signals of the first and second original A/D converter 117 into a digital image signal SI. The 

image signals may be generated, and the values of the digital image signal SI is stored on a storage medium 

sensitivities of the respective image signal components 65 118, such as a magnetic disk. Thereafter, the radiation 

of the first and second original image signals, which image stored on the second stimulable phosphor sheet 

represent picture elements in the corresponding radia- 104B is read out in the same manner as that described 

tion images, may be calculated on the basis of the un- above. The digital image signal S2 representing the 
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radiation image stored on the stimulable phosphor sheet S, a radiation image can be obtained in which noise has 
104B is stored on the storage medium 118. been minimized and which has good image quality. 

Thereafter, a superposition processing is carried out Alternatively, the sensitivities Ski and Sk2 of the two 
on the image signals Si and S2. FIG. 16 shows an appa- image signals SI and S2 may be detected for the respec- 
ratus for carrying out the superposition processing. 5 tive image signal components of the two image signals 
First, the image signals SI and S2 are read from an SI and S2 which represent image information stored at 
image file 118A and an image file 118B in the storage respective picture elements in the corresponding radia- 
medium 118 and fed into a superposition operating cir- tion images. How the sensitivity Ski of each of the 
cuit 119. The superposition operating circuit 119 calcu- image signal components of the image signal SI, which 
lates the sensitivity Ski of the image signal SI and the 10 represent image information stored at respective picture 
sensitivity Sk2 of the image signal S2. elements in the corresponding radiation image, is calcu- 

As indicated by line 70 shown in FIG. 11, the sensi- late d will oe described hereinbelow. 
tivity Ski and the sensitivity Sk2, which are calculated First > the sensitivity Ski and a latitude Gpl of the 
in this embodiment, are expressed in terms of the loga- unage signal SI are calculated from the image signal SI 
rithmic values of the sensitivity Ski' and the sensitivity 15 ^ method disclosed in U.S. Pat. No. 4,638,162. 
Sk2' detected from the image signals. Specifically, the For example, in cases where the sensitivity Ski is calcu- 
following formulas obtain. lated as being 2 and the latitude Gpl is calculated as 

being 2, the sensitivity Ski and the image signal SI have 
5Jti=iog skr, the relationship shown in FIG. 12. In cases where the 

20 picture element values are 10 bits, Smin==0, and 
s*2=Iog Ski' Smax= 1,023. When the x-y coordinates are set on the 

stimulable phosphor sheet 104A as shown in FIG. 13, 

After the sensitivity Ski and the sensitivity Sk2 are the sensitivity Skl(x,y) at an arbitrary point having the 
calculated, the image signals SI and S2 are weighted coordinates (x,y) on the stimulable phosphor sheet 
with the sensitivity Ski and the sensitivity Sk2 and are 25 104A can be calculated with the formula 
then added to each other. The weighted addition is 

represented by the formula _ Smax + Smin (21) 

S=(Sk2/(Skl+Sk2))XSlMSkl/{Skl+Skl))XS2 (10) Skl(x.y) - Smax - Smm X ~Gp~ + Skl 

30 

In this manner, an addition signal S is obtained from the wherein Sl(x,y) represents the value of the image signal 
two image signals SI and S2. Formula (10) corresponds component of the image signal SI, which image signal 
to Formula (9), in which Nl 2 is replaced by the sensitiv- component corresponds to an arbitrary point having the 
ity Skl, and N2 2 is replaced by the sensitivity Sk2. coordinates (x,y) on the stimulable phosphor sheet 
Therefore, noise contained in the addition signal S is 35 

minimized. 3to the manner described above, the sensitivity 

In cases where the weighted addition is carried out S^OwO corresponding to each of the picture elements 
by using the values of the sensitivities Skl' and Sk2' m tlie ra <^ at ion image, which was stored on the stimula- 
before being converted into the sensitivities Skl and ble P n ospkor sheet IMA, is calculated with Formula 
Sk2 having the logarithmic values, the addition signal S 40 ( 12 ^ Thereafter, the sensitivity Sk2(x,y) corresponding 
may be generated with the formula to of ^ P icture dements in the radiation image, 

which was stored on the stimulable phosphor sheet 
5= (log s&r/Oog SJfci'+iog Sft2O)x.si+0og 1MB, is calculated from the image signal S2. The addi- 

Skr/Qog SAi +iog skT))xsi (iff) tion signal S is then generated by carrying out the calcu- 

45 lations expressed as 

In cases where graduations are used for the sake of 
convenience for the sensitivity Sk shown in FIG. 11 S(x,y)=sk2(*yy(Ski(Ky)+sk2(x,y'))xsi+Ski(x- 
(e.g., in cases where graduations on the horizontal axis *y<^Cv)+SWfcj»)xfi <23) 

Sk" shown in FIG. 11 are used), the relationship be- T . « „ , , t 

tween Sk" and Sk' is expressed as 50 ^ the manner described above, the values of the 

sensitivities, from which the weight factors are deter- 
sjf a =4x io 4 -- 5 *" mined, axe calculated for the respective image signal 

components of the image signals, which represent pic- 
Therefore, in such cases, the addition signal S may be ture elements in the corresponding radiation images, 
generated with the formula 55 Therefore, noise can be imnimized for each of the pic- 

ture elements in the resulting image such that the image 
fc(4-5Jk2")/((4-5Jfcl")+(4-5Jt2"))x5H-(4-^ quality of the image may not be adversely affected by 

*i")/((4-s*i")+(4-S*2"))xS2 < 10 «) variations in the amount of radiation, which impinged 

upon the recording medium, such as the stimulable 
After the addition signal S is obtained in the manner 60 phosphor sheet, after being irradiated to the object, for 
described above, the addition signal S is fed into an different positions on the recording medium, 
image processing circuit 120 and subjected therein to Also, unsharp signals of the two image signals may be 
image processing, such as gradation processing and generated, and the values of the sensitivities of the re- 
frequency processing. The addition signal S, which has spective image signal components of the two image 
been obtained from the image processing circuit 120, is 65 signals, which represent picture elements in the cone- 
fed into an image reproducing apparatus 121 and used sponding radiation images, may be calculated on the 
during reproduction of a radiation image. When the basis of the unsharp signals of the two image signals. In 
radiation image is reproduced from the addition signal suc h cases, the values of the sensitivities of the respec- 
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tive image signal components of the two image signals, 
which represent picture elements in the corresponding 
radiation images, can be obtained such that the sensitivi- 
ties may not be adversely affected by minute changes 
and noise in the image signals. 5 

In the aforesaid embodiments of the superposition 
processing method for radiation images in accordance 
with the present invention, the weight factors to be used 
during the weighted addition are determined from the 
sensitivity Ski of the image signal SI and the sensitivity 10 
Sk2 of the image signal S2. Alternatively, noise Nl and 
noise N2 contained in the two image signals SI and S2 
may be detected, and the weighted addition may be 
carried out by using noise Nl and noise N2 in accor- 
dance with Formula (9). 15 

Also, in the aforesaid embodiments of the superposi- 
tion processing method for radiation images in accor- 
dance with the present invention, the sensitivity is cal- 
culated with the method disclosed in U.S. Pat. No. 
4,638,162 and Japanese Unexamined Patent Publication 20 
No. 61(1986>280163. Alternatively, any of other meth- 
ods, e.g., a method using a neural network, may be 
utilized for this purpose. 

Additionally, in the aforesaid embodiments of the 
superposition processing method for radiation images in 25 
accordance with the present invention, the superposi- 
tion processing is carried out on the image signals, 
which have been detected from the stimulable phosphor 
sheets. The superposition processing method for radia- 
tion images in accordance with the present invention is 30 
also applicable when the superposition processing is 
carried out on image signals obtained from other types 
of recording media, such as image signals detected from 
sheets of X-ray film or image signals obtained from an 
image intensifier. 35 

What is claimed is: 

1. A method for forming an energy subtraction image 
comprising the steps of: 

i) after two kinds of radiation images of an object are 
formed with two kinds of radiation having differ- 40 
ent energy levels, the object being constituted of a 
plurality of tissues exhibiting different levels of 
radiation absorptivity with respect to the two kinds 
of radiation having different energy levels, 

irradiating said two kinds of radiation images with 45 
image detection signals and detecting first and sec- 
ond original image signals based on said image 
detection signals and representing the two kinds of 
said radiation images, 

ii) generating a superposition image signal by carry- 50 
ing out a weighted addition expressed as 
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XSOl 



55 



wherein SOI represents said first original image signal, 
Nl represents noise contained in said first original 
image signal, S02 represents said second original image 
signal, and N2 represents noise contained in said second 
original image signal, 60 

iii) generating a first image signal, which represents a 
first image primarily composed of patterns of first 
tissues of said object, by subtracting said first origi- 
nal image signal and said second original image 
signal from each other, 65 

iv) generating a first smoothed image signal by 
smoothing said first image signal, said first 
smoothed image signal representing a first 



smoothed image in which noise components of said 
first image have been reduced or eliminated, and 
v) generating a second image signal by subtracting 
said superposition image signal and said first 
smoothed image signal from each other, said sec- 
ond image signal representing a second image pri- 
marily composed of patterns of second tissues of 
said object 

2. A method for forming an energy subtraction image 
as defined in claim 1 wherein said superposition image 
signal is generated by carrying out a weighted addition 
expressed as 

5=(5«/(5Jkl+5A2))xS01+(5«/(SJfcl+5«)>- 
XSOl 

wherein SOI represents said first original image signal, 
Ski represents a sensitivity of said first original image 
signal, S02 represents said second original image signal, 
and Sk2 represents a sensitivity of said second original 
image signal. 

3. A method for forming an energy subtraction image 
as defined in claim 1 wherein said superposition image 
signal is generated by calculating a sensitivity Ski of 
each of image signal components of said first original 
image signal, which represent picture elements in the 
corresponding radiation image, and a sensitivity Sk2 of 
each of image signal components of said second original 
image signal, which represent picture elements in the 
corresponding radiation image, on the basis of said first 
original image signal SOI and said second original 
image signal S02, and 

carrying out a weighted addition on the image signal 
components of said first original image signal and 
said second original image signal, which represent 
corresponding picture elements in the two kinds of 
said radiation images, said weighted addition being 
expressed as 

S~(Sk2/(Skl+Sk2))xSO%+(SkV(Skl+Sk2)y 
XSOZ 

4. A method for forming an energy subtraction image 
as defined in claim 1 wherein said superposition image 
signal is generated by generating an unsharp signal of 
said first original image signal SOI and an unsharp 
signal of said second original image signal S02, 

calculating a sensitivity Ski of each of image signal 
components of said first original image signal, 
which represent picture elements in the corre- 
sponding radiation image, and a sensitivity Sk2 of 
each of image signal components of said second 
original image signal, which represent picture ele- 
ments in the corresponding radiation image, on the 
basis of said unsharp signal of said first original 
image signal SOI and said unsharp signal of said 
second original image signal S02, and 

carrying out a weighted addition on the image signal 
components of said first original image signal and 
said second original image signal, which represent 
corresponding picture elements in the two kinds of 
said radiation images, said weighted addition being 
expressed as 

S-(SkV{Skl +SJB)) xSOX+(SkV(Skl+Sk2)y 
XS02 

5. A method for forming an energy subtraction image 
as defined in claim 1 wherein the two kinds of said 
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radiation images have been stored on stimulable phos- 
phor sheets. 

6. A method for forming an energy subtraction image 
as defined in claim 5 wherein each of said first and 
second original image signals is obtained by exposing 5 
each of said stimulable phosphor sheets to stimulating 
rays, which cause said stimulable phosphor sheet to 
emit light in proportion to the amount of energy stored 
thereon during its exposure to the radiation, and photo- 
electrically detecting the emitted light, 10 

7. A method for forming an energy subtraction image 
as defined in claim 6 wherein said stimulating rays are a 
laser beam. 

8. A method for forming an energy subtraction image 

as defined in claim 1 wherein the two kinds of said 15 
radiation images have been recorded on sheets of photo- 
graphic film. 

9. A method for forming an energy subtraction image 
comprising the steps of: 

i) after two kinds of radiation images of an object are 20 
formed with two kinds of radiation having differ- 
ent energy levels, the object being constituted of a 
plurality of tissues exhibiting different levels of 
radiation absorptivity with respect to the two kinds 

of radiation having different energy levels, 25 
irradiating said two kinds of radiation images with 
image detection signals and detecting first and sec- 
ond original image signals based on said image 
detection signals and representing the two kinds of 3Q 
said radiation images, 

ii) generating a superposition image signal by carry- 
ing out a weighted addition expressed as 
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wherein SOI represents said first original image signal, 
Nl represents noise contained in said first original 
image signal, S02 represents said second original image 43 
signal, and N2 represents noise contained in said second 
original image signal, 

iii) carrying out a first process for generating a first 
image signal, which represents a first image primar- 
ily composed of patterns of first tissues of said 45 
object, by subtracting said first original image sig- 
nal and said second original image signal from each 
other, 

iv) thereafter carrying out a second process, which 
comprises the steps of: 50 

a) generating a first smoothed image signal by 
smoothing said first image signal, said first 
smoothed image signal representing a first 
smoothed image in which noise components of 
said first image have been reduced, and 55 

b) generating a second image signal by subtracting 
said superposition image signal and said first 
smoothed image signal from each other, said 
second image signal representing a second image 
primarily composed of patterns of second tissues 60 
of said object, and 

v) thereafter carrying out a third process, which 
comprises the steps of: 

a) generating a second smoothed image signal by 
smoothing said second image signal, said second 65 
smoothed image signal representing a second 
smoothed image in which noise components of 
said second image have been reduced, and 



b) generating a new first image signal by subtract- 
ing said superposition image signal and said sec- 
ond smoothed image signal from each other, said 
new first image signal representing a new first 
image primarily composed of the patterns of said 
first tissues of said object 

10. A method for fonning an energy subtraction 
image as defined in claim 9 wherein said superposition 
image signal is generated by carrying out a weighted 
addition expressed as 

S={Skl/(Skl+Skl))xSOl+(Skl/(Skl+Sk2)y 
XS02 

wherein SOI represents said first original image signal, 
Ski represents a sensitivity of said first original image 
signal, S02 represents said second original image signal, 
and Sk2 represents a sensitivity of said second original 
image signal. 

11. A method for forming an energy subtraction 
image as defined in claim 9 wherein said superposition 
image signal is generated by calculating a sensitivity 
Ski of each of image signal components of said first 
original image signal, which represent picture elements 
in die corresponding radiation image, and a sensitivity 
Sk2 of each of image signal components of said second 
original image signal, which represent picture elements 
in the corresponding radiation image, on the basis of 
said first original image signal SOI and said second 
original image signal S02, and 

carrying out a weighted addition on the image signal 
components of said first original image signal and 
said second original image signal, which represent 
corresponding picture elements in the two kinds of 
said radiation images, said weighted addition being 
expressed as 

xsoi 

12. A method for forming an energy subtraction 
image as defined in claim 9 wherein said superposition 
image signal is generated by generating an unsharp 
signal of said first original image signal SOI and an 
unsharp signal of said second original image signal S02, 

calculating a sensitivity Ski of each of image signal 
components of said first original image signal, 
which represent picture elements in the corre- 
sponding radiation image, and a sensitivity Sk2 of 
each of image signal components of said second 
original image signal, which represent picture ele- 
ments in the corresponding radiation image, on the 
basis of said unsharp signal of said first original 
image signal SOI and said unsharp signal of said 
second original image signal S02, and 

carrying out a weighted addition on the image signal 
components of said first original image signal and 
said second original image signal, which represent 
corresponding picture elements in the two kinds of 
said radiation images, said weighted addition being 
expressed as 

S~(SJO/(Skl+Sk2)) xSOX+(Skl/<&a+Sk2)y 

xsoi 

13. A method for forming an energy subtraction 
image as defined in claim 9 wherein the two kinds of 
said radiation images have been stored on stimulable 
phosphor sheets. 
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14. A method for forming an energy subtraction 

image as defined in claim 13 wherein each of said first wherein SI represents the first image signal, Nl repre- 
and second original image signals is obtained by expos- sents noise contained in the first image signal, S2 repre- 
ing each of said stimulable phosphor sheets to stimulat- sents the second image signal, N2 represents noise con- 
ing rays, which cause said stimulable phosphor sheet to 5 tained in the second image signal, and S represents the 
emit light in proportion to the amount of energy stored addition signal. 

thereon during its exposure to the radiation, and photo- 23. A superposition processing apparatus for radia- 
electrically detecting the emitted light. tion images as defined in claim 22 wherein the radiation 

15. A method for forming an energy subtraction images have been stored on stimulable phosphor sheets, 
image as defined in claim 14 wherein said stimulating 10 24. A superposition processing apparatus for radia- 
rays are a laser beam. tion images as defined in claim 23 wherein each of the 

16. A method for forming an energy subtraction first and second image signals is obtained by exposing 
image as defined in claim 9 wherein the two kinds of each of said stimulable phosphor sheets to stimulating 
said radiation images have been recorded on sheets of rays, which cause said stimulable phosphor sheet to 
photographic film. 15 emit light in proportion to the amount of energy stored 

17. A superposition processing method for radiation thereon during its exposure to the radiation, and photo- 
images, comprising the steps of: electrically detecting the emitted light 

irradiating two kinds of radiation images of a single 25. A superposition processing apparatus for radia- 
object, each present on a medium, with image de- tion images as defined in claim 24 wherein said stimulat- 
tection signals and detecting, based on said image 20 ing rays are a laser beam. 

detection signals, first and second image signals, 26. A superposition processing apparatus for radia- 
each of which comprising a series of image signal tion images as defined in claim 22 wherein the radiation 
components and which represent corresponding images have been recorded on sheets of photographic 
picture elements in said radiation images; and film, 
adding the image signal components of the first and 25 27. A superposition processing method for radiation 
second image signals to each other in accordance images, comprising the steps of: 
with the following equation: irradiating two kinds of radiation images of a single 

object, each present on a medium, with image de- 
5=(A0 2 /M 2 +A2 1 ))x5i+(MVM 2 -f-^2 2 ))x52 tection signals and detecting, based on said image 

30 detection signals, first and second image signals, 
wherein SI represents the first image signal, Nl repre- each of which comprising a series of image signal 

sents noise contained in the first image signal, S2 repre- components and which represent corresponding 

sents the second image signal, N2 represents noise con- picture elements in said radiation images; and 

tained in the second image signal, and S represents the adding the image signal components of the first and 
addition signal. 35 second image signals to each other in accordance 

18. A superposition processing method for radiation with the following equation: 
images as defined in claim 17 wherein the radiation 

images have been stored on stimulable phosphor sheets. S~pk2/(Ski+sk2))xSt+(SkV(Ski+s*2))xSZ 

19. A superposition processing method for radiation 

images as defined in claim 18 wherein each of the first 40 wherein SI represents the first image signal, Ski repre- 
and second image signals is obtained by exposing each sents a sensitivity of the first image signal, S2 represents 
of said stimulable phosphor sheets to stimulating rays, the second image signal, Sk2 represents a sensitivity of 
which cause said stimulable phosphor sheet to emit light tnc second image signal, and S represents the addition 
in proportion to the amount of energy stored thereon signal. 

during its exposure to the radiation, and photoelectri- 45 28. A superposition processing method for radiation 
cally detecting the emitted light. images as defined in claim 27 wherein the radiation 

20. A superposition processing method for radiation images have been stored on stimulable phosphor sheets, 
images as defined in claim 19 wherein said stimulating 29. A superposition processing method for radiation 
rays are a laser beam. images as defined in claim 28 wherein each of the first 

21. A superposition processing method for radiation 50 311(1 second image signals is obtained by exposing each 
images as defined in claim 17 wherein the radiation of said stimulable phosphor sheets to stimulating rays, 
images have been recorded on sheets of photographic which cause said stimulable phosphor sheet to emit light 
film. in proportion to the amount of energy stored thereon 

22. A superposition processing apparatus for radia- during its exposure to the radiation, and photoelectri- 
tion images comprising: 35 cally detecting the emitted light. 

means for irradiating two kinds of radiation images of 30* A superposition processing method for radiation 
a single object, each present on a medium, with images as defined in claim 29 wherein said stimulating 
image detection signals and detecting, based on ravs are a laser beam. 

said image detection signals first and second image 31. A superposition processing method for radiation 
signals, each of which is made up of a series of & images as defined in claim 27 wherein the radiation 
image signal components and which represent cor- images have been recorded on sheets of photographic 
responding picture elements in said radiation im- film. 

ages; and 32. A superposition processing apparatus for radia- 

addition means for adding the image signal compo- ti° n images comprising: 
nents of the first and second image signals to each 65 means for irradiating two kinds of radiation images of 
other in accordance with the following equation: a single object, each present on a medium, with 

image detection signals and detecting, based on 
5=(JV2 2 /M J +A0 1 ))x5i-f <MViVl 2 +A2 3 ))xS2 said image detection signals first and second image 
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signals, each of which is made up of a series of 
image signal components and which represent cor- 
responding picture elements in said radiation im- 
ages; and 

addition means for adding the image signal compo 5 
nents of the first and second image signals to each 
other in accordance with the following equation: 

S=<Sk2/(Skl+Sk2))xSlHSkV(Skl+SkZ))xS2 

10 

wherein SI represents the first image signal, Ski repre- 
sents a sensitivity of the first image signal, S2 represents 
the second image signal, Sk2 represents a sensitivity of 
the second image signal, and S represents the addition 
signal. 

33. A superposition processing apparatus for radia- 
tion images as defined in claim 32 wherein the radiation 
images have been stored on stimulable phosphor sheets. 

34. A superposition processing apparatus for radia- 
tion images as defined in claim 33 wherein each of the 
first and second image signals is obtained by exposing 
each of said stimulable phosphor sheets to stimulating 
rays, which cause said stimulable phosphor sheet to 
emit light in proportion to the amount of energy stored 25 
thereon during its exposure to the radiation, and photo- 
electrically detecting the emitted light. 

35. A superposition processing apparatus for radia- 
tion images as defined in claim 34 wherein said stimulat- 
ing rays are a laser beam. 30 

36. A superposition processing apparatus for radia- 
tion images as defined in claim 32 wherein the radiation 
images have been recorded on sheets of photographic 
film. 

37. A superposition processing method for radiation 35 
images, comprising the steps of: 

i) irradiating two kinds of radiation images of a single 
object, each present on a medium, with image de- 
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40. A superposition processing method for radiation 
images as defined in claim 39 wherein said stimulating 
rays are a laser beam. 

41. A superposition processing method for radiation 
images as defined in claim 37 wherein the radiation 
images have been recorded on sheets of photographic 
film. 

42. A superposition processing apparatus for radia- 
tion images comprising: 

means for irradiating two kinds of radiation images of 
a single object, each present on a medium, with 
image detection signals and detecting, based on 
said image detection signals, first and second image 
signals, each of which comprising a series of image 
signal components and which represent corre- 
sponding picture elements in said radiation images; 
and 

means for generating an addition signal based on the 
image signal components of the first and second 
image signals, comprising: 

i) an operation device for calculating a sensitivity Ski 
of each of the image signal components of the first 
image signal and a sensitivity Sk2 of each of the 
image signal components of the second image sig- 
nal on the basis of the first image signal SI and the 
second image signal S2; and 

ii) an addition device for adding the image signal 
components of the first image signal and the second 
image signal in accordance with the following 
equation to obtain an addition signal S; 

S={Sk2/(Skl+Sk2))*SlHSkV(Skl+Sk2))XS2. 

43. A superposition processing apparatus for radia- 
tion images as defined in claim 42 wherein the radiation 
images have been stored on stimulable phosphor sheets. 

44. A superposition processing apparatus for radia- 
tion images as defined in claim 43 wherein each of the 



2512^ ° n 40 fct ™« second * obtained by exposing 

detection signals, first and second image signals, 40 each of said stimulable phosphor 



each of which comprising a series of image signal 
components and which represent corresponding 
picture elements in said radiation images; 
ii) calculating a sensitivity Ski of each of the image 
signal components of the first image signal and a 
sensitivity Sk2 of each of the image signal compo- 
nents of the second image signal on the basis of the 
first image signal SI and the second image signal 
S2; and 

in) adding the image signal components of the first 
image signal and the second image signal in accor- 
dance with the following equation to obtain an 
addition signal S: 

S=(Sk2/(Skl+Sk2))XSl+(SkV{Skl+Sk2y)XS2. 
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38. A superposition processing method for radiation 
images as defined in claim 37 wherein the radiation 
images have been stored on stimulable phosphor sheets, $q 

39. A superposition processing method for radiation 
images as defined in claim 38 wherein each of the first 
and second image signals is obtained by exposing each 
of said stimulable phosphor sheets to stimulating rays, 
which cause said stimulable phosphor sheet to emit light 65 
in proportion to the amount of energy stored thereon 
during its exposure to the radiation, and photoelectri- 
cally detecting the emitted light 



sheets to stimulating 
rays, which cause said stimulable phosphor sheet to 
emit light in proportion to the amount of energy stored 
thereon during its exposure to the radiation, and photo- 
electrically detecting the emitted light 

45. A superposition processing apparatus for radia- 
tion images as defined in claim 44 wherein said stimulat- 
ing rays are a laser beam. 

46. A superposition processing apparatus for radia- 
tion images as defined in claim 42 wherein the radiation 
images have been recorded on sheets of photographic 
film. 

47. A superposition processing method for radiation 
images, comprising the steps of: 

i) irradiating two kinds of radiation images of a single 
object, each present on a medium, with image de- 
tection signals and detecting, based on said image 
detection signals, first and second image signals, 
each of which comprising a series of image signal 
components and which represent corresponding 
picture elements in said radiation images; 

ii) generating an unsharp signal of the first image 
signal SI and an unsharp signal of the second image 
signal S2; 

iii) calculating a sensitivity Ski of each of the image 
signal components of the first image signal and a 
sensitivity Sk2 of each of the image signal compo- 
nents of the second image signal on the basis of said 
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unsharp signal of the first image signal SI and said 
unsharp signal of the second image signal S2; and 
iv) adding the image signal components of the first 
image signal and the second image signal in accor- 
dance with the following equation to obtain an 
addition signal S: 

S=(Sk2/(Skl+Sk2))XSl+{Skl/{Skl+SkX))XS2. 



10 



48. A superposition processing method for radiation 
images as defined in claim 47 wherein the radiation 
images have been stored on stimulable phosphor sheets. 

49. A superposition processing method for radiation 
images as defined in claim 47 wherein each of the first 15 
and second image signals is obtained by exposing each 
of said stimulable phosphor sheets to stimulating rays, 
which cause said stimulable phosphor sheet to emit light 

in proportion to the amount of energy stored thereon 
during its exposure to the radiation, and photoelectri- 20 
cally detecting the emitted light 

50. A superposition processing method for radiation 
images as defined in claim 49 wherein said stimulating 
rays are a laser beam. 

51. A superposition processing method for radiation 2 * 
images as defined in claim 47 wherein the radiation 
images have been recorded on sheets of photographic 
film. 

52. A superposition processing apparatus for radia- 
tion images comprising: 

means for irradiating two kinds of radiation images of 
a single object, each present on a medium, with 
signals image detection and detecting, based on 
said image detection signals, first and second image 35 
signals, each of which comprising a series of image 
signal components and which represent corre- 
sponding picture elements in said radiation images; 
and 
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means for generating an addition signal based on the 
image signal components of the first and second 
image signals, comprising: 

i) an operation device for generating an unsharp 
signal of the first image signal SI and an unsharp 
signal of the second image signal S2, and calcu- 
lating a sensitivity Ski of each of the image sig- 
nal components of the first image signal and a 
sensitivity Sk2 of each of the image signal com- 
ponents of the second image signal on the basis 
of said unsharp signal of the first image signal SI 
and said unsharp signal of the second image sig- 
nal S2; and 

ii) an addition device for adding the image signal 
components of the first image signal and the 
second image signal in accordance with the fol- 
lowing equation to obtain an addition signal S: 

S=(Sk2/{Skl+Sk2))xSl+(SkV(Skl+SkX))xSl 

53. A superposition processing apparatus for radia- 
tion images as defined in claim 52 wherein the radiation 
images have been stored on stimulable phosphor sheets. 

54. A superposition processing apparatus for radia- 
tion images as defined in claim 53 wherein each of the 
first and second image signals is obtained by exposing 
each of said stimulable phosphor sheets to stimulating 
rays, which cause said stimulable phosphor sheet to 
emit light in proportion to the amount of energy stored 
thereon during its exposure to the radiation, and photo- 
electrically detecting the emitted light 

55. A superposition processing apparatus for radia- 
tion images as defined in claim 54 wherein said stimulat- 
ing rays are a laser beam. 

56. A superposition processing apparatus for radia- 
tion images as defined in claim 52 wherein the radiation 
images have been recorded on sheets of photographic 
film. 

***** 
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